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Systematic deviations from standard Fermi-liquid behavior have been widely observed and documented in several 
classes of strongly correlated metals. For many of these systems, mounting evidence is emerging that the anoma- 
lous behavior is most likely triggered by the interplay of quenched disorder and strong electronic correlations. In 
this review, we present a broad overview of such disorder-driven non-Fermi-liquid behavior, and discuss various 
examples where the anomalies have been studied in detail. We describe both their phenomenological aspects as 
observed in experiment, and the current theoretical scenarios that attempt to unravel their microscopic origin. 
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I. INTRODUCTION 

Soon after the discovery of high temperature superconductivity in the cuprates it was reaUzed that the optimally doped com- 
pounds showed clear deviations from Fermi Uquid behavior. This was not the first time that a real system seemed to violate 
Landau's paradigm. The heavy fermion superconductor UBeia had previously been found to show anomalous normal state 
behavior (Cox, 1987). Since then, the quest for such systems has been vigorous and numerous puzzling compounds have been 
discovered. However, a coherent theoretical framework with which to understand this behavior has not been obtained. One 
of the few situations where there seems to be a general agreement as to the inadequacy of Fermi liquid theory is the critical 
region governed by a zero temperature quantum phase transition. Here, deviations from Landau's predictions seem to be the 
norm, although a general theory of this behavior is still lacking. However, the quantum critical non-Fermi liquid (NFL) requires 
fine tuning to a special location in the phase diagram. In many of the known compounds, NFL behavior seems to be better 
characterized as the property of a whole phase instead of just a phase boundary. 

Interestingly, NFL behavior in disordered systems rests on a much firmer theoretical basis. Even an exactly solvable model is 
known where diverging thermodynamical properties can be calculated ( McCoy, 1969; McCoy a ndJWu, 1968) . In addition, this 
anomalous behavior persists over an entire region of the phase diagram and is not confined to a specific point. These so-called 
Griffiths-McCoy phases had been originally proposed in disordered classical systems, where the accompanying singularities are 
rather weak ( Griffiths, 1969). The quantum version of this phenomenon is, however, much stronger and has now been established 
theoretically in many models. Furthermore, their properties are not believed to be a peculiarity of the one-dimensional exactly 
solvable model. Their relevance to real compounds has also been advocated in many systems, from doped semi-conductors to 
disordered heavy fermion compounds. The emergence of NFL behavior in a disordered context is thus seen to be fairly natural. 

The question of NFL behavior in various clean situations has been extensively discussed in the literature. It is the purpose of 
this review to attempt to bring together in a single article a number of theoretical models and analyses in which disorder-induced 
NFL behavior has arisen. Although reference will be made to the most important experimental findings, the focus will be on 
theory rather than on experiments. For the latter, we will refer the reader to other reviews. 

The issue of disorder-induced NFL behavior was discussed for the first t ime in connection with the anomalous the rmody- 
namic responses near the metal-insulator transition in doped semiconductors dPaalanen and Bhat3.ll99lLISarachiklll995l) . In this 
context, the random singlet phase of Bhatt and Lee was an important pioneering work ('Bhat t and LeeLll98lLll982r Since then, 
anomalous behavior has been observed in various disordered heavy fermion compounds, where local moments are well formed 
even in the clean limit (Stewart, 2001) . More recently, tw o-dimensional electron systems in semiconductor heterostructures 
have also shown intriguing behavior ( Abrah ams et a/.Lf2001 ). Finally, intrinsic and extrinsic heterogeneities in high temperature 
superconductors have led to the consideration of these disorder effects in that context as well. These systems, therefore, have 
become the main focus of the theoretical efforts and will take center stage in this review. 

We will attempt to emphasize the general aspects of the various contexts in which NFL behavior arises. As mentioned above, 
Griffiths phases play a central role in many of these models. Although they come in different guises, all Griffiths phases seem 
to share a few common generic features, which we will describe. Once these features are present, most physical quantities 
can be immediately obtained. The question of distinguishing among the different microscopic mechanisms is therefore rather 
subtle. However, a few criteria can be posed and perhaps used to make this distinction. We hope these will serve as a guide to 
experimentalists. 

This review is organized as follows. In Section |II] we briefly summarize the main ideas of Fermi liquid theory as applied 
to both clean and disordered systems. Section |ni| is devoted to a general discussion of the phenomenology of the main non- 
Fermi liquid systems. Subsection IIII. Al focuses on NFL behavior of clean systems, whereas subsections IIII.BI IlII.CI and llll.DI 
describe the anomalous behavior of doped semiconductors, disordered heavy fermion systems, and the glassy regime close to 
the metal-insulator transitions. The main theoretical approaches are reviewed in Section HVl Subsection IIV. Al is devoted to 
results on the disordered Hubbard and Kondo/Anderson lattice models. We highlight the physics of local moment formation 
in a disordered environment dlV.A.H . the random singlet phase ( IIV.A.2I I. the Kondo disorder model ( IIV.A.3l and IIV.A.41 . the 
electronic Griffiths phase close to the disorder-induced metal-insulator transition ( IIV.A.5I IIV.A.6I and IIV.A.7> . and the inco- 
herent transport near the two-dimensional metal-insulator transition ( IIV.A.8> . The general topic of magnetic Griffiths phases 
is discussed in subsection llV.BI We first review their well-established properties in the case of insulating magnets, where we 
emphasize the important role played by the infinite randomness fixed point ( IIV.B.U . The generic mechanism behind magnetic 
Griffiths phases is described in subsection II V.B .21 the important effect of dissipation on the dynamics of the Griffiths droplets 
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is considered in subsection lIV.B.3l and further remarks on the appHcability of this scenario to real metalHc systems are made 
in subsection II V.B.4I Section HV CI is devoted to the physics of glassy dynamics and ordering in metals. We wrap up with 
conclusions and open questions in SectionlVl 



II. GENERAL PREDICTIONS OF FERMI LIQUID THEORY AND ITS LIMITATIONS 
A. Fermi liquid theory for clean systems 

Interacting homogeneous systems of fe rmions in three dim ensions are believed to be described at low energies and long wave 
lengths by Landau's Fermi liquid theory llLandaulfl957albl) . Although there is no general rigorous iustificatio n for it. Landa u 
was able to use diagrammatic arguments to show that his general framework is at least internally consistent (Landau, '19591). 
More recently, these early demonstrations were put on a firmer basis in the modern framework of the renorma lization group, in 
which the asymptotic correctness of the theory can be shown llBenfatto and Gallavottii[l99(ilShankailll994l) . In addition, our 
confidence in its correctness is boosted by simple model systems where controlled calculations can be carried out, such as a 
dilute Fermi gas with short-range interactions ( Abrikosov and Khalatnikov, 1958; Galitskii, 1958), the electron-phonon system 
jMigdalU 19581) . and the very dense degenerate electron plasma dAbrikosovl 1 1 962) . More importantly, its direct experimental 
verification in the prototypical Fermi liquid system, the normal state of liquid ^He between the onset of super fluidity at 3 
mK and about 100 mK, can be taken as one of the most striking exam ples of Landau' s theory at work (see, e. g., ( lLeggettl|l975h ). 
Finally, Fermi liquid theory, properly generalized to charged systems ( ISihnlfT958albh . forms the foundation of our understanding 
of the behavior of electrons in metals. 

There are excellent reviews of Fermi liq uid theory to which the reader is referred to jAbrikosov et a/.Lri975tlBavm and Pethicld 
11991; Leegett, 1975; Pi nes and Noziere i 11965^. We will content ourselves with a brief outline of its main assumptions and 
predictions. Fermi liquid theory starts by recognizing that the low energy excitations of a Fermi sea, say an added electron at a 
wave vector k with — k > kp and k — kp <C kp, where /cf is the Fermi wave vector (we assume a rotationally invariant 
system, for simplicity), have a very long lifetime. This is due to restrictions imposed by energy and momentum conservation 
and the blocking of further occupation of the Fermi sea by the Pauli exclusion principle. In this case, the lifetime can be 
shown to be Tk ~ (fc — kp)^ ■ The stability of these quasi-particle excitations allows us to ignore their decay at first and treat 
them (perturbatively) only at a later stage. The second ingredient of the theory is sometimes called "adiabatic continuity": 
the excitations of the interacting fermions are in one-to-one correspondence to those of a non-interacting Fermi gas, in other 
words, the slow "switching-on" of the interactions in a Fermi gas does not alter the nature of the excitation spectrum. These two 
assumptions led Landau to write the total energy of weakly excited states of the system as a functional of the various occupation 
numbers of states labeled by momentum k and spin projection a ~ ±, the quasi-particle states. More precisely, he used the 
deviations (5nk,cr from their occupations in the ground state (in this review, we use units such that ft = 1 and ks = 1) 

k,(T 

+ ^ /(k, cr;k',cr')(5nk.<T<5nk'.<T'- (1) 

k,(T;k' ,cr' 

Here, Eq is the (usually unknown) ground state energy, e (k) ^ vp {k — kp) is, the quasi-particle dispersion, parametrized by 
vp — kp/m* , where to* is the effective mass of the quasi-particle, usually distinct from its ba re, free-electron v alue. Note that 
the Fermi wave vector kp \?, not modified by interactions, a fact known as Luttinger's theorem (Luttinger, 1960). The last term, 
quadratic in the occupations, incorporates the self-consistent interactions among the quasi-particles and is easily seen to be of 
the same order as the second term for weakly excited states. 

Since only momenta close to the Fermi surface are needed for low-energy excitations, we can take, to leading order, |k| = 
|k'| = kp in the last term of Eq. Q and then / only depends on the angle 9 between k and k'. Moreover, in the absence of an 
applied magnetic field, spin rotation invariance imposes additional constraints 

/(k,a;k',a') = .r(0) + (Ta'r (0). (2) 

Finally, it is convenient to decompose the angular dependence in Legendre polynomials 

vpf^- [0) = ^ (2; + 1) F^'^'Pi (cos 9) , (3) 

where vp = Vm*kp /tt^ is the total density of states at the Fermi level and V is the volume of the system. The dimensionless 
constants Fj"'" are known as Landau coefficients. 
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We thus see that the low-energy sector of the interacting system is completely parametrized by the effective mass to* and the 
Landau coefficients, in terms of which important physical properties can be written. For example, the low temperature specific 
heat is linear in temperature with a coefficient determined by the effective mass 

^ = '-m*k,. (4) 
T 3 

The magnetic susceptibility is a constant at low temperatures and related to the Landau parameter Fq 

X(T)^4;^, (5) 
TT^ (1 + Fg) 

where /is is the Bohr magneton. The compressibility k ^ —V {dP/ dV), which can be accessed through the sound velocity, is 
related to Fq 

Ott^to* 

(6) 



4(1 + ^0)' 

There are several other predictions which come out of the basic framework outlined above, especially with regard to collective 
excitations (zero sound, plasmons, spin waves) and transport properties, but we will not dwell on them, referring the interested 
reader to the available reviews. We will only note that the quadratic dependence on its excitation energy of the quasi-particle 
lifetime leads to a quadratic temperature dependence of the resistivity at the lowest temperatures. Moreover, in the T — > limit 
the resistivity tends to a constant value determined by lattice imperfections and extrinsic impurities, such that the Fermi liquid 
prediction is 

p{T) = po+AT\ (7) 

The thermal conductivity of a Fermi liquid is proportional Ta (T) (the Wiedemann-Franz law), where a (T) — 1/p (T) is the 
electric conductivity. 

The application of Fermi liquid theory to metals helps explain why a strongly interacting electron gas in a clean crystal can 
behave almost like a free electron gas and conduct heat and electricity so well at low temperatures. Interestingly, very strong 
electronic correlations can put the theory to test at the most extreme circumstances. In particular, in a certain class of compounds 
containing rare-earth or actinide elements such as Ce or U (to be reviewed later in Section lTlI.A.H . the effective mass is observed 
to be 2 to 3 orders of magnitude greater than the electron mass. These compounds were therefore named "heavy fermions". 
Nevertheless, even in these cases, Fermi liquid theory provides a valid description in a number of cases, although a growing 
number of exceptions have been discovered and investigated in recent years. Furthermore, the idea of adiabatic continuation (the 
one-to-one correspondence between the excitations of a simpler reference system and another one of interest) has proven fruitful 
in contexts which go far beyond Landau's original proposal. In particular, this philosophy has been extended to the superfluid 
phases of '^He (.LeggetL.lQTS.) and to nuclear physics (.Migdal. .19,67). More importantly for the subject of this review, Fermi 
liquid ideas form the basis of much of our understanding of interacting electrons in disordered metals, as we will expand upon 
in the next Section. 



B. Fermi liquid theory for disordered systems 

While the original formulation and much of the later work on Fermi liquid theory concentrated on clean metals, the relevant 
physical principles have a much more general validity. This framework is flexible enough to be also applicable not only in 
presence of arbitrary forms of randomness, but even for finite size systems such as quantum dots, mol ecules, atoms, or atomic 
nuclei. In electronic systems, the first systematic studies of the interplay of interactions and disorder dLee and Ram akrishnanl 
Ii9§5) 

emerged only in the last 25 years or so. Most progress was achieved in the regime of weak disorder, where controUed 
many-body calculations are possible using the disorder strength as a small parameter 



1. Drude theory 

To lowest order in the disorder strength, one obtains the semi-classical predictions of the Drude theory 
(iLee and Ramakrishnanl fT985h . where the conductivity takes the form 

a « do = , (8) 
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where n is the carrier concentration, e the electron charge and m its band mass. According to Matthiessen's rule, the transport 
scattering rate takes additive contributions from different scattering channels, viz. 

rtr' = ^eT' + r-e' (T) + t'^' (T) + • • • . (9) 

Here, tJ^ is the elastic scattering rate (describing impurity scattering), and t^J-{T), t~^{T),..., describe inelastic scattering 
processes from electrons, phonons, etc. It is important to note that in this picture the resistivity p = cr^^ is a monotonically 
increasing function of temperature 

p{T) « + AT", (10) 

since inelastic scattering is assumed to only increase at higher temperatures. The residual resistivity po = (j{T — 0)~^ is 
thus viewed as a measure of impurity (elastic) scattering. Within this formulation, thermodynamic quantities are generally not 
expected to acquire any singular or non-analytic corrections due to impurity scattering. 

Drude theory encapsulates a very simple physical picture. It implicitly assumes that the itinerant carriers undergo many 
collisions with unspecified scattering centers, but that these scattering events remain independent and uncorrected, justifying 
Matthiessen's rule. At this level, inelastic scattering processes are therefore assumed to be independent of impurity scattering, 
and thus assume the same form as in standard Fermi liquid theory, e.g. t~i^ ^ T^, etc. This simplifying assumption is better 
justified at higher temperatures, where inelastic scattering processes erase the phase memory of the electrons, and suppress 
quantum interference processes arising from multiple scattering events. 



2. Perturbative quantum corrections 

At weak disorder, systematic corrections to the Drude theory were found jLee and RamakrishnanlflQSSh to consist of several 
additive terms, 

Cr = (To + <5(7tu; + (^CTint, (11) 

corresponding to the so-called "weak localization" and "interaction" corrections. These "hydrodynamic" corrections are domi- 
nated by infrared singularities, i.e., they acquire non-analytic contributions from small momenta or equivalently large distances, 
and which generally lead to an instability of the paramagnetic Fermi liquid in two dimensions. Specifically, the weak localization 
corrections take the form 



i-(d-2) _ T-{d-2) 



(12) 



where I — vftk the mean free path, d is the dimension of the system, and Lrh is the length scale over which the wave functions 
are coherent. This effective system size is generally assumed to be a function of temperature of the form Lth ^ Tp/"^, where the 
exponent p depends on the dominant source of decoherence through inelastic scattering. The situation is simpler in the presence 
of a weak magnetic field where the weak localization corrections are suppressed and t he lea ding dependence comes from the 
interaction corrections first discovered by Altshuler and Aronov ( Altshuler and Aro nov.ll979h 

Sa.^t = ^(ci - C2F.)(Tt)(^-2)/2. (13) 
n 

Here, Ci and C2 are constants, and is an interaction amplitude. In d ~ 3, this leads to a square-root singularity 5<Tint ~ "s/T, 
and in d = 2 to a logarithmic divergence Saint ^ In(Tr). These corrections are generally more singular than the temperature 
dependence of the Drude term, and thus they are easily identified experimentally at the lowest temperatures. Indeed, the T^/^ 
law is commonly observed (Lee and Ramakrishnan , 1985) in transport experiments in many disordered metals at the lowest 
temperatures, typically below 500 mK. 

Similar corrections have been predicted for other physical quantities, such as the tunneling density of states and, more impor- 
tantly, for thermodynamic response functions. As in Drude theory, these quantities are not expected to be appreciably affected 
by noninteracting localization processes, but singular contributions are predicted from interaction corrections. In particular, cor- 
rections to both the spin susceptibility and the specific heat coefficient 7 — Cy / T were expected to take the general forms 
iJt^ ^ ^ X{d.-2)/2^ ^jjj j.jjjjg {hree dimensions 



X-Xo + m^VT, (14) 
7 = 7o + m^VT, (15) 
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where and are constants. 



As in conventional Fermi liquid theories, these corrections emerged already when th e interactio ns were treated at the lowest, 
Hartree-Fock level, as done in the approach of Altshuler and Aronov ( Altshuler and Ar onoyl [19791) . Higher order corrections in 
the interaction amplitude were first incorporated by Finkelshtein (Finkel'stein, 1983, 1984), demonstrating that the predictions 
remained essentially unaltered, at least within the regime of weak disorder. In this sense, Fermi li quid theory has been generalize d 
to weakly disordered metals, where its predictions have been confirmed in numerous materials llLee and Ramakrishnanlfl985l) . 



3. Scaling theories of disordered Fermi liquids 

In several systems, experimental tests were extended beyond the regime of weak disorder, where at le ast at face value, the per- 
turbative predictions seem of questionable relevance. Interestingly, a number of transport experiments ( iLee and Ramakrishnant 
[1985: Paalanen andBhatt, 1991; Sarachik, 1995) seemed to indicate that some predictions, such as the T^/^ conductivity law, 
appear to persist beyond the regime of weak disorder At stronger disorder, the system approaches a dis order-driven metal- 
insulator transition. Since the ground-breaking experiment of Paalanen, Rosenbaum, and Thomas in 1980 llRosenbaum et al.{ 
[1980), it became clear that this is a continuous (second order) phase transition ( Paalanen ef a/. , 1982), which bears many 
similaritie s to conventional critical phenomena. This important observation has sparked a veritable avalanche of exper- 
imental ( Lee and Ramakrishnan', '1985', 'Paalanen an d Bhatl Il99lt LSarachiki 11995*) and theoretical tAbraham s et al.[ '1979^ 
ISchaffer and Wesner, 1980; Wegner, 1976, 1979) works, most of which have borrowed idea s from studies of second order 
phase transitions. Indeed, many experimental results wer e interpreted u sing scaling concepts jLee and Ramakrishnanl 11983) . 
culminating with the famous scaling theory of localization jAbrahams et aL.1979.) . 

The essential idea of these approaches focuses on the fact that a weak, logarithmic instability of the clean Fermi liquid arises 
in two dimensions, suggesting that d ~ 2 corresponds to the lower critical dimension of the problem. In conventional critical 
phenomena, such logarithmic corrections at the lower critical dimension typically emerge due to long wavelength fluctuations as- 
sociated wit h spontaneously broken continuous symmetry. Indeed, early work of Wegner ISchaffer and Wegneil[l98(ilWegneil 
Il974ll979 h emphasized the analogy between the localization transition and the critical behavior of Heisenberg magnets. It 
mapped the problem onto a field theoretical nonlinear cr-model and identified the hydrodynamic modes leading to singular cor- 
rections in d = 2. Since the ordered (metallic) phase is only marginally unstable in two dimensions, the critical behavior in 
d > 2 can be investigated by expanding around two dimensions. Technically, this is facilitated by the fact that in dimension 
d = 2+e the critical value of disorder W for the metal-insulator transition is very small (Wc ^ e), and thus can be accessed using 
perturbative renormalization group (RG) approa ches in direct an alogy to the procedures developed for Heisenberg magnets. 

In this approach ( Abrahams et al. .1979: .Schaffer and Wegner. .1980.: .Wegner. .1976. ■1979) . conductance is identified as the 
fundamental scaling variable associated with the critical point, which is an unstable fixed point of the RG flows. In this picture, 
temperature scaling is obtained from examining the system at increasingly longer length scales Lxh ~ T^/"^, which follows 
from the precise form of the RG flows. In the metallic phase, under rescaling the conductance g oo, (corresponding to the 
reduction of effective disorder), and the long distance behavior of all correlation functions is controlled by the approach to the 
stable fixed point at = 0. In other words, the leading low temperature behavior of all quantities should be identical to that 
calculated at infinitesimal disorder This scaling argument therefore provides strong justification for using the weak disorder 
predictions throughout the entire metallic phase, provided that the temperature is low enough. 



4. Fermi liquid near localization transitions 

In the following years, these ideas were extended with a great deal of effort in t he formulation of a scal i ng theory of interacting 
disordered electrons by Finkelshtein ('Finkel'stein', '1983") and many followers jBehtz and Kirkpatricld 1 19941 ICastellani et al\. 
^984). While initially clad in an apparent veil of quantum field theory jargon, these theories were later given a simple physical 
interpretation in terms of Fermi liquid ideas (Castella ni et a/.L Il987; Kotliar, 1987) for disordered electr ons. Technical details of 
these theories a re of co nsiderable complexity and the interested reader is referred to the original literature ( iBelitz and KirkpatrickL 
[1994; Castella ni et al\. [l984; Finkel'stein, 1983, 1984). Here we just summarize the principal results, in order to clarify the 
constraints imposed by these Fermi liquid approaches. 

Within the Fermi liquid theory for disordered systems dCastellani et al l [T987[lKotliail[T987l) . the low energy (low temperature) 
behavior of the system is characterized by a small number of effective parameters, which include the diffusion constant D, the 
frequency renormalization factor Z, and the interaction amplitudes 7^ and 74. These quantities can also be related to the 
corresponding quasi-particle parameters which include the quasi-particle density of states 

PQ = Zpo, (16) 

and the quasi-particle diffusion constant 



Dq = DIZ ^ D/pQ. 



(17) 
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Here, po is the "bare" density of states which describes the noninteracting electrons. In the absence of intera ctions, the sing le- 
particle density of states is only weakly modified by disorder and remains noncritical (finite) at the transition ( IWegneiill98ll) . 

Using these parameters, we can now express the thermodynamic response functions as follows. We can write the compress- 
ibility 

Xc = ^=Pq[1-27«], (18) 

the spin susceptibility 

Xs = ^ilpQ[l-2-^t], (19) 

and the specific heat 

Cv = 2tt^ pQT/i. (20) 
In addition, we can use the same parameters to express transport properties such as the conductivity 

dn 

o = j-^D, = pqDq, (21) 
as well as the density-density and spin-spin correlation functions 



dpD(.q —lOJ DsQ'' — ■"' 



Here, we have expressed these properties in terms of the spin and charge diffusion constants, which are defined as 

_ D _ D 

^^^Z(l-27.)' ^^=Z(l-270" ^^^^ 

Note that the the quantity D is not the charge diffusion constant Dc that enters the Einstein relation, Eq. J21> . As we can 
write a = poD, and since po is not critical at any type of transition, the quantity D (also called the "renormalized diffusion 
constant") has a critical behavior identical to that of the conductivity cr. We also mention that the quasi-particle diffusion constant 
Dq = D/Z has been physically interpreted as the heat diffusion constant. 



5. Critical scaling 

Finally, we discuss the scaling behavior of observables in the crit ical region. In particular, if the scaling description 
llBehtz and KirknatricldllQQllCastellani J^ll984lFinkersteinHl983[ll984 ) is valid, the conductivity can be written as 

a{t, T) = 6-(''-2) (b^/'^ t, b'T). (24) 

Here, b is the length rescaling factor, and t = {n ~ nc)/nc is the dimensionless distance from the transition. We have also 
introduce d the c orrelation length exponent ly and the "dynamical exponent" z. This expression, first proposed by Wegner 
([Abrahams et a/.L fT979: Schaffer and Wegner, 1980; Wegner, 1976, 1979), is expected to hold for "regular"types of transitions, 
where the critical values of the interaction amplitudes remain finite. The conductivity exponent p in a {T = Q) ~ can be 
obtained by working at low temperatures and choosing b ~ t^'^ . We immediately see that 

(7(T) ^t^0,(T/r^), (25) 

where <j)cr{x) = f„{\,x) and 

p^{d-2)u, (26) 
a relation known as "Wegner scaling". Finite temperature corrections in the metallic phase are obtained by expanding 

<^„(a;) w l + ax", (27) 

giving the low temperature conductivity of the form 



a{t,T)K<Jo{t)+m„{t)T". 



(28) 



8 



Here, adt) ~ f^, and merit) ~ ^A'-ai^z Since the scaling function 4>a{x) is independent of t, the exponent a must take 
an universal value in the entire metallic phase, and thus it can be calculated at weak disorder, giving a — 1/2. This scaling 
argument provides a formal justification for using the predictions from perturbative quantum corrections as giving the leading 
low temperature dependence in the entire metallic phase. Note, however, that according to this result the prefactor ma{t) is not 
correctly predicted by perturbative calculations, since it undergoes Fermi liquid renormalizations which can acquire a singular 
form in the critical region near the metal-insulator transition. 

The temperature dependence at the critical point (in the critical region) is obtained if we put t ~ Q (n — Uc), and choose 
6 = T-i/^ We get 

cr(i = 0,T) ~ t(''"2)A^ (29) 

An analogous argument can be used llCastellani and Castro[ll98^ for the specific heat coefficient 7 = Cv/T [it is important 
not to confuse this quantity with the interaction amplitude 74] 

7(t, T) = (6^/'^ t, b'^T). (30) 

Choosing h = t^" and expanding in T we find 

7(i,T)«7„(i) + m^(i)ri/2, (31) 

with 7o(t) ~ t-'^, m^{t) - f-('=+''^/2). Even though 7^ = 7(t, T = 0) can become singular at the transition, it is expected to 
remain finite away from the transition it 7^ 0). W e stress that the specific heat exponent k and the dynamical exponent z are not 
independent quantities. In fact, Castellani et al. dCastellani and Castrcilll98^ have proved that within Fermi liquid theory, the 
following relation is obeyed 

z = d+-. (32) 

V 

Similar conclusions apply also to other quantities, such as the spin susceptibility Xs, which should also remain bounded away 
from the critical point and acquire universal T^l"^ corrections within the metallic phase. Finally, the compressibility Xc = dn/d^ 
is generically expected to remain nonsingular (finite) at the transition , since n{^) is expected to be a smooth function, except at 
special filling fractions. For example, if one approaches ( lMotllfT990l) a band or a Mott insulator, the compressibility vanishes as 
a precursor of the gap opening at the Fermi surface. 

All the above expressions are quite general, and can be considered to be a phenomenological description dCastellani et al\. 
Il987^ of disordered Fermi liquids. On the other hand, these relations tell us nothing about the specific values of the Landau 
parameters, or how they behave in the vicinity of the metal-insulator transition. Perturbative renormalization group calculations 
jBelitz and Kirkpatrick, 1994; Castellani et al., 1984; Finkel'stein, 1983, 1984) based on the 2 + e expansion have been used to 
make explicit predictions for the values of the critical exponents and scaling functions for different universality classes. Despite 
a great deal of effort invested in such calculations, the predictions of these perturbative RG approaches have not met almost 
any success in explaining the experimental data in the critical region of the metal-insulator transition. We should emphasize, 
though, that limitations associated with these weak-coupling theories do not invalidate the potential applicability of Fermi liquid 
ideas per se. On the other hand, even in their most general form, these Fermi liquid considerations predict that thermodynamic 
response functions such as x 7 remain finite at T = even within the disordered metallic phase. All experiments that find 
a more singular behavior of these quantities should therefore be regarded as violating the Fermi liquid theory - even when it is 
properly generalized to disordered systems. 



C. Typical departures from Fermi liquid theory In real systems 

Fermi liquid theory describes the leading low energy excitations in a system of fermions. These, as we have seen, are weakly 
interacting quasiparticles characterized by several Landau parameters. Most remarkably, its validity is by no means limited to 
systems with weak interactions. In several materials, most notably heavy fermion systems, these many body renormalizations 
are surprisingly large (e.g. m*/m ^ 1000), yet Fermi liquid theory is known to apply at the lowest temperatures. 

On the other hand, we should make it clear that precisely in such strongly correlated systems the temperature range where 
Fermi liquid theory applies is often quite limited. The "coherence temperature" T* below which it applies is typically inversely 
proportional to the effective mass enhancement, and thus can be much smaller then the Fermi temperature. Above T* all 
physical quantities are dominated by large incoherent electron-electron scattering, and Fermi liquid theory simply ceases to be 
valid. In many heavy Fermion systems T* is found to be comparable to the single-ion Kondo temperature, thus to typically be 
of order 10- WK. 

Another interesting class of materials, where clear departures from Fermi liquid theory have been observed at sufficiently 
high temperatures, are systems close to the Mott transition. These include the transition metal oxides such as {yi-xCrx)20^ 
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Figure 1 Temperature dependence of the resist ivity at d i fferen t pressures for a two-dimensional organic charge transfer k — (BEDT — 
TTF)2Cu[N{CN)2]Cl, following Ref. iLime lette e? a/.Ll2003h . The data (circles) are compared to the DMFT predictions (diamonds), with 
a pressure dependence of the bandwidth as indicated. 

and chalcogenides such as NiS2^xSex (for a review see Ref. dlmada et all \l99§ti ). In very recent work, similar behavior 
was observed in anothe r Mott system, a two-dimensional organic charge transfer salt k — {BEDT — TT F)2Cu[N {C N)2]Cl 
dLimelette et a/.LE003l) . In this material the system can be pressure-tuned across the Mott transition, and on the metallic side the 
resistivity follows the conventional law below a crossover temperature T* ~ bi)K. Above this temperature transport crosses 
over to an insulating-like form, reflecting the destruction of heavy quasiparticles by strong inelastic scatter ing (see Fig. [Hi . 
Clearly, Fermi liquid theory does not apply above this coherence scale, but alternative DMFT approaches dLimelette et al\. 
|2003) proved very successful in quantitatively fitting the data over the entire temperature range. 

The incoherent metallic behavior could be even more important in disordered systems, where the coherence temperature 
T* can be viewed as a random, position-dependent quantity r*(x). If this random quantity is characterized by a sufficiently 
broad probability distribution function P{T*), then Fermi liquid theory may not be applicable in any temperature regime. 
In the following, we examine several physical systems where these phenomena are of potential importance, and then describe 
theoretical efforts to address these issues and produce a consistent physical picture of a disorder-driven non-Fermi liquid metallic 
state. 



III. PHENOMENOLOGY OF NON-FERMI LIQUID BEHAVIOR AND EXPERIMENTAL REALIZATIONS 
A. NFL behavior in correlated systems with weak or no disorder 

1 . Heavy fermion materials near magnetic quantum critical points 

Heavy fermion systems have been the focus of intense study over the last three decades as prime examples of strongly 
correlated electronic behavior. These are metallic systems with an aiTay of localized moments formed in the incompletely-filled 
/-shells of rare-earth or actinide elements. An important requirement for heavy fermion behavior is that the /-shell electrons 
should be sufficiently close to a valence instability and therefore should have a fairly low ionization energy in the metallic host. 
This restricts the interesting behavior to a few elements, most often Ce and U, but also Yb, Pr, Sm and some other less common 
cases. The proximity to a valence instability promotes the enhancement of the hybridization between the /-shell electrons and 
the conduction bands. This, in turn, leads to an antiferromagnetic (AFM) interaction between the /-shell moment and the local 
conduction electron spin density through essentially a super-exchange mechanism. This interaction is otherwise much too weak 
in the lanthanide or actinide series. The tendency towards local singlet formation through the Kondo effect serves to strongly 
suppress the ubiquitous magnetic order seen in most other intermetallics with localized /-moments. The hybridization of the 
metallic caiTiers with the virtual bound states in the strongly correlated /-shell leads to large effective mass renormalization 
factors (of order lO'^ — 10'^), whence the name "heavy fermions". There are many excellent reviews of heavy fermion physics to 
which the reader is referred to (Grewe and Steglich, 1991; Hewson, 1993; Lee et ai, 1986; Stewart, 1984). 

A great part of the interesting physics of these compounds is a result of competing tendencies: localization versus delocal- 
ization of the /-electrons and magnetic order versus Kondo singlet formation. This complex interplay provides the background 
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c Aui showing the logarithmic divergence at criticality. Data are from Ref. iBogenbereer and v. Lohnevseri 



where a wide variety of low temperature phases arise: there are metals, insulators ( lAeppli and Fisld.ri992ft . and superconductors 
( Heffner and Norman, 1996). Furthermore, many exotic phenomena are also found, such as small moment antiferromagnetism 
( Buyers. .1996) . coexistence of superconductivity and magnetism, and unconventional superconductivity with more than one 
phase jHeffaer and Norman, 1996). 

Despite a bewildering zoo of correlated behavior, in recent years a common theme has been the focus of attention. In many 
systems, the Neel temperature is low enough that it can be tuned to zero by an external parameter, such as pressure, chemical 
pressure (by alloying with an element with a different ionic radius) or applied magnetic field. The zero temperature phase 
transition between an antiferromagnet and a paramagnet as a function of the external parameter is called a quantum phase 
transition. If this happens to be a second order phase transition, the system will exhibit a diverging correlation length and 
we expect the critical behavior to be clas sified in "universality classes", much like thermal second order phase transitions 
(IContinentinol 119941 1200 it |Sachdevl Il999h . Indeed, many heavy fermion systems can be tuned in just this way and a general 
theory of this so-called quantum critical point (QCP) has been sought vigorously. The possibility of such a quantum phase 
transition in heavy fermion systems is generally attr ibuted to Doniach, who discussed it early on in the context of a one- 
dimensional effective model ( lDoniachlfl977t) (see also dContinentino et aD.ll98 9)). The associated phase diagram is thus usually 
referred to as the Doniach phase diagram. It is a natural consequence of the above mentioned competition between Kondo singlet 
formation and AFM order. 

Besides the natural interest in a catalogue of the possible universality classes of such quantum phase transitions, the low 
temperature region in the vicinity of the zero temperature quantum critical point is observed to be characterized by strong 
deviations from Landau's Fermi liquid theory. A recent exhaustive review of the experimental data on many heavy fermion 
system s showing NFL behavior, including but not restricted to those where QCP physics seems relevant, can be found in cStewarl 
l200ll) . In many cases, the following properties are observed: 

• A diverging specific heat coefficient, often logarithmically, C (T) /T ~ 70 log {Tq/T) (see Fig.|2}. 

• An anomalous temperature dependence of the resistivity, p (T) ^ po + AT", where a < 2. 

• An anomalous Curie-Weiss law, x"^ (T) ^ Xo^ + CT^, where /3 < 1. 



Despite an intensive theoretical onslaught, a coherent theoretical picture is still lacking dColeman et all 1200 it iMillisl 
[n)99). The natural theoretical description of the metallic QCP would seem to be one where the critical modes are s ubject 
to the interactions generated in the presence of a Fermi sea, as in the old paramagnon theory ("Berk and Schrieffer, 'l966j 
[opniach and Enaelsbera, 1966; Izuyar na et aL , 1963) . In particular, the low-frequency dependence of the RPA susceptibility 
leads to the overdamping of the critical modes, which can decay into particle-hole pairs (Landau damping). As emphasized 
by Hertz, this increases the effective dimensionality of the critical theory by the dynamical exponent z: def f = d + z ( Hert J, 
|T976) (see also ( Beal-Monod and Maki, 1975)). The dynamical exponent z determines the relation between energy and length 
scales E ^ L^^ . Landau damping leads Xo oj ^ —ik^iz = 3) in clean itinerant ferromagnets and lo ^ —iq^ (z — 2) in 
clean itinerant antiferromagnets (q being the deviation from the ordering vector), the latter being the case of interest in heavy 
fermion materials. As a result, three-dimensional systems are generically above their upper critical dimension, the dimension 
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above which mean field exponents become exact (equal to 4 in the usual thermal case for Ising and Heisenberg symmetries). In 
the effective critical theo ry, the critical modes are th e more weakl y interacting th e longer the length scales, rendering the theory 
asymptotically tractable dContinentinol 1 1 994t iHertzL fl 976 : MiUisl l 1 993t iMorviaL ITqSS ). This is the weak coupling spin density 
wave approach, sometimes named the Hertz-Millis scenario. In particular, at the AFM QCP the specific heat is non-singular, 
C (T) /T ~ 7o — a^/T. Furthermore, scattering off the incipient AFM order is singular only along lines on the Fermi surface 
which are connected by the AFM ordering wave vector Q ("hot spots"). These are effectively short-circuited by the remainder 
of the Fermi surface ("cold spots"), which are hardly affected by criticality, leading in clean samples to a Fermi liquid response 
p (T) ~ po + (Hlubina and Rice, 1995). These results cannot explain the observed properties of heavy fermion systems 
close to quantum criticality. Another consequence of the higher effective dimensionality is the fact that the theory does not obey 
hyperscaling. Hyperscaling is usually associated with the Josephson scaling law: 2 — a = vd, where a is the specific heat critical 
exponent and v is the correlation length exponent (Goldenfeld, 1992). In a quantum phase transition, d should be replaced by 
the effective dimensionality de//- However, the most important consequence of the violation of hyperscaling is the absence of 
uj/T scaling in dynamical responses which, in the AFM case, should involve the combination lo/T'^/'^ instead (Coleman et ali 
[200 1). However, uj/T scaling is observed in a quantum critical heavy fermion system CeCu5.9Auo.i (|Schroder ef a/., 199^ 
I2OOOI) . These results present major difficulties for the description of the NFL behavior observed near a QCP in heavy fermion 
systems. 

There have been proposals for amending the Hertz-Millis scenario in a n effort to account fo r the experimental results. In par- 
ticular, we mention the effect of disorder close to a quantum critical point ( Rosch, 1999, 2000) and the possible two-dimensional 
character of the spin fluctuation spectrum (Mathur et ai, 1998; Roschgf ai, 1997; Schr oder et ai, 2000; Stockert et al , 199 8j). 
Disorder acts by relaxing the momentum conservation which restricts strong scattering to the "hot lines", effectively enhancing 
their influence on transport. We stress that in this case disorder plays only an ancillary role and is not the driving mechanism for 
NFL behavior On the other hand, a two-dimensional spin fluctuation spectrum places the system at the upper critical dimension 
and is able to account for some of the observed critical exponents. None of these approaches, however, is able to account for 
all the available data. Alternatively, radical departures from the spin density wave critical theory have been proposed. In one 
of these, the long time dynamics of the locali zed moment acquires a non-trivial power-law dependence, whereas the spatial 
fluctuations retain the usual mean-field form (ISi et a/.Ll200ll) . This has been dubbed 'local quantum criticality' and is a natural 
scenario to incorporate the lu/T scaling with non-trivial exponents observed over much of the Brillouin zone. Its current form, 
however, still relies on a two-dimensional spin fluctuation spectrum. In another approach, there is a phase transition from the 
usual paramagnetic heavy Fermi liquid state to a non-trivial state where the localized spins form a 'spin liquid' weakly coupled 
to the conduction sea ( Senthil ef a/. , 2003, 2004). An interesting consequence of the latter proposals would be a drastic red uction 
of the Fermi surface volume across the transition, presumably measurable through the Hall constant (IColeman et a/.llToOll) . This 
is still an open arena where perhaps new ideas and new experiments will be necessary before further progress can be achieved. 

2. Marginal Fermi liquid behavior of high Tc superconductors 

High temperature superconductors (HTS) have taken the center stage of modern condensed matter theory ever since their 
discovery in 1987. Most of their features continue to defy theoretical understanding, and even the origin of superconducting 
pairing in these materials remains controversial. Most remarkably, many features of the superconducting state seem to be less 
exotic than the extremely unusual behavior observed in the normal phase. In fact, HTS materials presented one of the first well 
documented examples featuring deviations from Fermi liquid phenomenology, initiating much theoretical activity and debate. 
On a phenomenological level, behavior close to "optimal doping" (the regime where the superconducting transition temperature 
Tc is the highest) seems to display "marginal Fermi liquid" (MFL) behavior The most prominent feature observed in experiments 
is the linear resistivity (see Fig.|3}, which at optimal dop ing persists in an enormous temperature range from a few Kelvin all 
the way to much above room temperature (Takaai et al, 1992). A phenomenological model describing the MFL behavior of 
cuprates has been put forward a long time ago ( Varma et al... .1989 ), but its microscopic origin remains highly controversial. 
To our knowledge no microscopic theory has so far been able to provide a microscopic explanation for this puzzling behavior, 
despite years of effort and hundreds of papers published on the subject. 

Much theoretical and experimental effort over the years concentrated on unraveling the microscopic origin of the supercon- 
ducting pairing, a mechanism that, one hopes, would also explain the puzzling features of the normal state. Most proposed 
scenarios concentrated on exotic many-body mechanisms in an assumed homogeneous conductor, thus ignoring the role of dis- 
order or the possible emergence of some form of random ordering. Many excellent reviews exist detailing the current status of 
these research efforts and, given our focus on the effects of disorder, will not be elaborated here. In the following, we follow 
Ref. ( Panagop oulos and Dobrosavlievicll2005l) . and briefly review those experimental and theoretical efforts suggesting that the 
homogeneous picture may be incomplete. 

The first evidence of random ordering in HTS was obtained from the observation of low temperature spin-glass ordering in 
the pseudo-gap phase. At doping concentration x larger than a few percent, antiferromagnetic ordering is suppressed, but the 
short range magnetic order persists. The low-field susceptibility displays a cusp at temperature T ^ Tg and a thermal hysteresis 
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Figure 3 Normal state in-plane resistivity of optimally doped cuprate La2-xSrxCu04 (circles). The striking linear temperature dependence 
is observed around optimal doping, persis ting far above room temperature. Deviations from this behavior are seen both below (squares) and 
above optimal doping. Data are from Ref. iTakagi ef fl/.lll992h . 
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Figure 4 Phase diagram of the archetypal HTS, following Ref. I'Panagopoulos and Dobros avlievid .f200.'j^. Tn is the Neel temperature, Tf 
and Tg the onset of short range freezing to an electronic glass, and Tc the superconducting transition temperature. At a; < Xsc the material is 
a glassy insulator. At Xsc < x < Xopt a microscopically inhomogeneous conducting glassy state emerges, with intercalated superconducting 
and magnetic regions. At a; = Xopt the system experiences a quantum glass transition and at a; > Xopt the material transforms into a 
homogeneous metal with BCS-like superconducting properties. The superfluid density is maximum at a:: = Xopt- The crossover scales T* and 
Tm characterizing normal-state transport (see text for details), vanish at the quantum glass transition. 



below it, characteristic of a spin glass transition. At T < Tg the mat erial displays memory effects like "traditional" spin glasses 
and is described by an Edwards-Anderson order parameter dChou et ai, 199^. Such magnetic measurements are not possible 
in the superconducting regime {x > Xsc), but muon spin relaxation (/iSR) has been successful in identifying the freezing of 
electronic mom ents under the superconducting dome of various HTS systems (Kaniael et ai, 2002; Panagopoulos et ai. 20031 
l2002tlSanna et ai. 2004). From these studies one may conclude that in this regime glassiness coexists with superconductivity 
on a microscopic scale throughout the bulk of the material. Most remarkably, the observed spin-glass phase seems to end at a 
quantum critical point precisely at optimal doping (x = Xopt), suggesting that non-Fermi liquid behavior in the normal phase 
may be related to the emergence of glassy ordering in the ground state (Fig.|3. 

The coiTelation between spin order and charge transport is further emphasized by experiments in high magnetic fields 
feoebinger et ai. 1996) where bulk superconductivity was suppressed, revealing information about low-T charge transport in 
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the normal phase. The resistivity data show a crossover in pab {T) from metallic to insulating-like behavior (resistivity minimum) 
at a characteristic temperature T* which, similarly to Tg, decreases upon doping, and seems to vanish at the putative quantum 
critical point. In addition, a crossover temperature T,„ alx > Xopt se parating marginal F ermi liquid transport at T > T,„ from 
more conventional metallic behavior at T < also seems to drop ( Naaib et al, 2003) to very small values around optimal 
doping (see Fig.|3i. At x > Xopt the gr ound state becomes metallic and homogeneous, with no evidence for glass iness or other 
form of nano-scale heterogeneity (Bala kirev et aZ.L l2()03l:lBoebinger et aZlllQQdlPanagopoulos et aZ.Ll2003[E00^ . Remarkable 
independent evidence that a QCP is found precisely at a; = Xopt is provided by the observation of a sharp change in the super- 
fluid density ns{0) ^ l/^ah(0) (where AQb(O) is the absolute value of the in-plane penetration depth). At x > Xopt, ns{0) is 
main ly doping independent (F ig.Bt, while the T-dependence is in good agreement with the BCS weak-coupling d-wave predic- 
tion ( Panagopoulos et ai\\2003t) . At dopings below the quantum glass transition ns{0) is rapidly suppressed (note the enhanced 
depletion near a; = 1 /8 precisely where T„ and T * are enhanced) and there is a marked departure of (T) from the canonical 
weak coupling curve (iPanagopoulos et oZT 12003). All these results provide strong evidence for a sharp change in ground state 
properties at a; = Xopt, and the emergence of vanishing temperature scales as this point is approached - just as one expects at a 
QCP. 

The formation of an inhomogeneous state below optimal doping is consistent with those theoretical scenarios that predict 
phase separation (.Dag otto, 2002: Gor'kov and Sokol, 1987: Kiyelson et al.,. 2003) at low doping. Coulomb interactions, how- 
ever, enforce charge neutrality and prevent jKivelsonef 0^120031) global phase separation; instead, the carriers are expected 
(ISchmalian and Wolvnest '2000*) to segregate into nano-scale domains - to form a stripe/cluster glass ('Schmalian and Wo lvnesl 
I2000h . As quantum fluctuations increase upon doping (Dobrosavlievic ef a/., 2003; Pastor and Dobrosavlievic, 1999), such a 
glassy phase should be eventually suppressed at a quantum critical point, around x = Xopt- These ideas find striking support 
in very interesting STM studies revealing nano-scale domains forming in the underdoped phase, as recently reported by several 
research groups. 

This scenario is consistent with another mysteriou s aspect of norm al state transport in the weakly underdoped regime 
(x ^ Xopt)- Here, DC transport has a much weaker l lBoebinger et ail 11996 ) (although still insulati ng-like) tem perature de- 
pendence. However, the observed logT resistivity upturn in this region has been shown (Boebinger to be incon- 
sistent with conventional localization/interaction corrections which could indicate an insulating ground state. Instead, estimates 
feeloborodov et al... 2003 ) reveal this behavior to be consistent with that expec ted for metallic droplet charging/tunnelling pro- 
cesses, as seen in quantum dots and granular metals ( iBeloborodov et all l2003l) . These results suggest that in this regime HTS 
are inhomogeneous metals, where conducting droplets connect throughout the sample, and a metal-insulator transition in the 
normal phase happens exactly at x — Xgc- At lower densities the conducting droplets remain isolated, and the material may be 
viewed as an insulating cluster or stripe glass. As carrier concentration increases they connect and the carriers are free to move 
throughout the sample, forming filaments or "rivers". This is, in fact, the point where free carriers emerge in Hall effect data 
iBalakirev et al., 2003) and phase coherent bulk superconductivity arises at a; > Xgc- This observation suggests that it is the 
inhomogeneous nature of the underdoped glassy region which controls and limits the extent of the superconducting phase at low 
doping. 

At this time we still do not know if the formation of such inhomogeneous states is indeed a fundamental property of HTS 
materials, or merely a by-product of strong frustration and extrinsic disorder Nevertheless, the emerging evidence seems com- 
pelling enough by itself, as it opens the possibility that inhomogeneities and glassy ordering may not be disregarded as the 
possible origin of non-Fermi liquid behavior in the cuprates. 

B. Disorder-driven NFL behavior in correlated systems: doped semiconductors 

Doped semiconductors llShklovskii and Efrosill98l have been studied for a long time, not the least because of their enormous 
technological applications, but also because of their relatively simple chemical composition allowing the possibility for simple 
theoretical modelling ( Shklovskii and Efros. 1984) . Typically, they are used to fabricate transistors and other devices, which 
essentially can be used as electrical switches in logical integrated circuits. For this reason, much attention has been devoted to 
understanding their behavior close to the metal-insulator transition (MIT) (Mott, 1990), which typically takes place when the 
Bohr radius Oo of the donor atoms becomes comparable to the typical carrier separation. Most applications are based on doped 
silicon devices, where the critical concentration n = ng ^ lO^^cm" ^. Since the Fermi energy of the carriers in this regime is 
fairly low {Tp k, IQQK) dPaalanen and Bhat?. '1991', 'Sarachik", '1995V the behavior at room temperatures is easy to understand 
using conventional solid-states theories (Paalanen and Bhatt, 1991; Sarachik, 1995), and as such it has been qualitatively and 
even quantitatively understood for many years. The evolution of the physical properties as a function of n, however, is smooth 
at elevated temperatures due to thermal activation. 
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Figure 5 Critical behavior of tlie conductivity extrapolated to T ^ for uncompensated Si:P. Note a sharp power-law behavior with exponent 
/i ~ 1/2, extending over a surprisingly large concentration range. The results reported are taken from the classic stress-tuned experiment of 
Paalanen, Rosenbaum, and Thomas liRosenbaum et aL.198Q,) . 



1. Metal-insulator transition 

A sharp metal-insulator transition is seen only at the lowest accessible temperatures, and it is in this regime where most of the 
interesting physics associated with quantum effects comes into play. 

In the last 30 years, a great deal of effort has been devoted to the study of detailed properties (Paalanen and Bhatt, '19911 



IS^achiL 1995J of doped semiconductors in the critical region near the MIT. Early experiments ( Paalanen and Bhatt^ 199j 
1 1199.'^ concentrated on transport measurements, and reported a temperature dependence of the conductivity of the form 



a{T) ^cTo + mVT, (33) 



where (Tq and m are parameters that generally depended on the carrier concentration n and the magnetic field B. Using these 
observations, one typically extrapolated these results to T = 0, and examined the critical behavior as the transition is ap proached. 
In agreement with scaling predictions ( Abrahams et al. Vl979^. 'Behtz and Kirkpatrick'. '1994^ 'Fin keFsteinl [l983l Il984) . power- 
law behavior was reported ( Paalanen and Bhatt, 1991; Rosenbaum et ai, 1980: Sarachik, 1995) of the form (see Fig.|5} 

a{T ^ 0) - (n - UcY. (34) 

The conductivity exponent was generally reported to take the value /i k, 0.5 for uncompensated materials (e.g. only donor 
dopants present as in Si:P or only acceptor s as in Si:B) in zero magnetic field. In the presence of a magnetic field experiments 
dPaalanen and Bhatllll99lllSarachikl ll995^ showed /I w 1, and the same behavior was reported in compensated materials (both 
donors and acceptors present, such as Si:P,B). These findings had the general form expected from Fermi liquid considerations 
and scaling approaches, althou gh the specific values for the conductivity exponent ^ were no t easy to understand from the 
available microscopic theories dBelitz and KirkpatrickL [l99j. It was, however, generally felt CLee and Ramakrishnanl Il985h 
that the observed compensation dependence reflects the enhanced role of the electronic correlations in the uncompensated case, 
where the dopant impurity band is half-filled, and one would expect a Mott metal-insulator transition in the absence of disorder. 



2. Thermodynamic anomalies 

If the proximity to the Mott transition (Mott, 1990) is an important consideratio n, then magnetism may be expected to have 
unusual features, since the electrons turn into local magnetic moments ('Mott'.ll990| ) in the Mott insulator. Indeed, many strongly 
correlated metals near the Mott transition (e.g. transition metal oxides ( Motg. [l990l) such as V2O3) s how large spin suscepti- 
bility and specific heat e nhancements, which are known to be dominated by local magnetic moment dOuirt and Markol ll97lL 
lUeand Maekawalll97ll) physics. 
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Figure 6 Temperature dependence of the normalized s uscep tibility for three different Si:P,B samples, with electron densities n/ric — 
0.58, 1.1, 1.8, as reported by Hirsch et al. iHirsch et ^^93). The low temperature behavior looks qualitatively the same in a very broad 
concentration interval, persisting well into the metallic side of the transition, in contrast to Fermi liquid predictions. 



To test these ideas and also to explore the vahdity of generalized Fermi liquid considerations dCastellani et al\. Il987h . sub- 
sequent experiments t urned to exploring the thermodynamic behavior of the system. On the insulating side, one expects 
jShklovskii and EfrosL flQSl the donor electrons to be tightly bound to ionic centers, thus behaving as spin 1/2 local mag- 
netic moments. Because of the overlap of the donor ele ctron wave-functions, antiferromagnetic exchange is generated between 
pairs of such magnetic moments, which takes the form llShklovskii and Efroslll984l) 

J(i?) Joexp{-i?/a}. (35) 

The parameters Jo and a (the Bohr radius) can be calculated with precision for all the known shallow impurity centers. Since the 
donor atoms are randomly distributed throughout the host matrix, the insulator can therefore be described as a random quantum 
antiferromagnet, which at low temperature can be expected to exhibit random spin freezing, i.e. spin glass behavior. Despite 
considerable effort, the search for such spin glass ordering in the insulating phase has proven unsuccessful ( Ouirt and Mark^ 
(197 it lUe and Maekawal 1197 ll) . Instead, a very unusual thermodynamic response has been observed, with a diverging spin 
susceptibility and sub-linear specific heat as a function of temperature, viz. 

x(r)^T"-\ (36) 

l{T) = ^ ^ r"-i. (37) 

The exponent a was found dPaalanen and Bhatl ll99lL ISarachikl 1 19951) to weakly depend on the details of specific system, 
generally taking the value a « 0.4 for uncompensated and a « 0.3 for compensated systems. Soon after the initial observations, 
this behavior on the insulating side was qualitatively and even quantitatively explained by the "random singlet" theory of Bhatt 
and Lee cBhatt and Lee, 1981, 1982) (see Section llVXIl . 



3. Two-fluid model 



The real surprise, however, came from the obser vation that essentially the same low temperature anomalies persisted 
dPaalanen et al\. 179881 ISchlager and Lohnevsenl Il997^ well into the metallic side of the transition. The data were fitted by a 
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phenomenological "two-fluid" model (" Paalanen et al[\\92,'A\Onin and Markd.ll97lh . which assumed coexistence of local mag- 
netic moments and itinerant electrons, so that (see Fig.|S} 

7/70 = m*/ml + [To/Tr-\ (38) 
X/Xo = m*/ml + (i{TjTr-\ (39) 

In these expressions, the first term is ascribed to itinerant electrons which are assumed to form a Fermi liquid with band mass m* 
and effective mass to*. The second term describes the contribution from localized magnetic moments, wh ere the deviation from 
the C urie law (a > 1) is viewed as a result of spin-spin interactions, similarly as in the B hatt-Lee theory llBhatt and Leeifl98]l 
Il982h . The exponent a was found jPaalanen et al\. Il988l ISchlaeer and Lohnevsenl Il99 7) to have a weak density dependence, 
and essentially take the same value as on the insulating side. The fitted values of the parameter To can be used to estimate 
the relative fraction of the electrons that form the local moments. According to these estimates, between 10% and 25% of the 
electrons contribute to the formation of these localized moments in the vicinity of the transition. 

These experimental findings represented drastic violations of the Fermi liquid predictions. Within the Fermi hquid picture 
(ICastellani et a/.L[l987l) . in any metal the local magnetic moments should be screened by conduction electrons through the Kondo 
effect or spin-spin interactions, so that 7 and x should remain finite as T ^ 0. Their precise value should be a function of the 
corresponding Fermi liquid parameters m*and (see Eqs. (0} and (|5j), and as such should have a singular behavior only as a 
phase transition is approached, not within the metallic phase. In contrast, the singularities observed in these materials display no 
observable anomaly (Hirsch et al , 1992) as one crosses from the metallic to the insulating side. In fact, the thermodynamic data 
on their own have such a weak density dependence, that just on their basis one could not even determine the critical concentration 
ric- This behavior is in sharp contrast to that observed in conventional Fermi liquids (e.g. clean heavy fermion compounds), 
where sharp anomalies in both the transport and the thermodynamic properties are clearly seen and well documented. 

Since the physics associat ed with precursors of l ocal m oment magnetism is typically associated with strong correlation ef- 
fects, it seems very likely tLee and Ramakrishnanl Il985l) that the surprising behavior of doped semiconductors has the same 
origin. If this is true, then the conventional weak-coupling approaches (Belitz and Kirkpatrick, 1994) to electronic correlations 
seem ill suited to describe this puzzling behavior These ideas seem in agreement with the original physical picture proposed 
by Mott (.Mott ,1990.) . where the Coulomb repulsion is viewed as a major driving force behind the metal-insulator transition. 
Indee d, even early attempts to a ddress the role of correlation effects based on Hubbard models by Milovanovic, Sachdev, and 
Bhatt ('Mi lovanovic et all . 1989) were sufficient to indicate the emergence of disorder-induced local moment formation around 
the transition region. This theory w as too simple to address more subtle questions such as the interac tion of such local mo- 
ments and the conduction electrons jBhgttand Fishe it 1 19921: iDobrosavlievic ef al. . Il992t iLakner et al.i 0L994), but later work 
jDobrosavlievic an d G.Kotlia'3.ll993lll994lTl997tiDobrosavhevic and Kotliailll998l) based on extended DMFT approaches pro- 
vided further support to the strong correlation picture. 



4. Asymptotic critical beliavior or mean-field scaling? 

One more important aspect of the experimental data is worth emphasizing. Several well defined features of these materials 
have been clearly identified, both with respect to the thermodynamic and the transport properties. These include a sharply 
defined power-law behavior of the conductivity, with an exponent fi being a strong function of the magnetic field or the deviation 
from half filling, and a smooth but singular thermodynamic response. All of these features are clearly defined over a very broad 
range of parameters covering, for example, in excess of 100% of the reduced concentration 5n = {n — nc)/nc. Such behavior 
provides a strong indication that the puzzhng observations are not associated with complications arising within the asymptotic 
critical region of a second-order ph ase transition. In contrast, the behavior in the immediate vicinity of the critical point is 
still a subject of some con troversy (Stupp et'all 1 1 9931 IWaffenschmidt et'ol. and may be dominated by extrinsic (e.g. 

self-compensation) effects d ltoh et a/.n2004liK 

Therefore, theories designed to capture only long wavelength and low energy excitations do not seem likely to be sufficient in 
describing the most puzzling features of the experiments. To our knowledge, none of the microscopic predictions of such "field 
theoretical" approaches ( Belitz an d Kirkpatrick,, 1994) have been experimentally verified; most observations seem drastically 
different from what one expects based on these theories. Such robust and systematic behavior observed in a broad parameter 
range is more reminiscent of mean-field behavior, which typically describes most of the experimental data around phase transi- 
tions, except in extremely narrow critical regions. Typically, one has to approach the critical point within a reduced temperature 
t = {T — Tc)/Tc of the order of 10^"^ or less in order to observe the true asymptotic critical behavior. Unfortunately, an appro- 
priate mean-field description is not available yet for metal-insulator transitions, despite years of effort. In our opinion, what is 
needed to correctly describe the observed non-Fermi liquid features are approaches that can explicitly address the interplay of 
strong correlations and disorder. 
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C. Disorder-driven NFL behavior in correlated systems: disordered heavy fermion systems 

Anomalous NFL behavior has been observed in a large number of non-stoichiom etric heavy ferrn ion compounds. An ample 
review of the experimental situation until 2001 can be found in Stewart's review (Stewart, '2001"). Since the main focus of 
this review is on the theoretical approaches to disordered NFL systems, we will not dwell much on the available data on these 
systems. We will, however, give a bird's eye view of the subject. 

We would first like to distinguish between the main systems of interest here and those whose detailed description requires 
the consideration of disorder effects but in which disorder is most likely not the driving mechanism of NFL behavior This is 
especially important in those cases where the proximity to a quantum critical point seems to be the origin of the anomalous 
behavior Indeed, the external tuning parameter in these systems is often chemical pressure, namely, the substitution of small 
amounts of an isoelectronic element with a slightly larger ionic radius with the aim of expanding the lattice. Non-stoichiometry 
is then unavoidable but its effects are considered secondary when compared with, say, the variation of the couplings between 
neighboring atoms. Some studies of these effects, especially on transport have been undertaken ( Rosch. 1999. .2000 ). A caveat 
is in order, however A well-known criterion due to Harris, establishes that if i^d < 2, where v is the correlation length critical 
exponent of the clean system, then disorder is a relevant perturbation, i. e., the actual critical behavior is modified in the dirty 
system cl I arris. 1974). The clean Hertz-Millis theory discussed in Section lTll.A. ll has a mean field exponent = 1/2 and the 
Harris criterion indicates that disorder is relevant. Therefore, even in this case it is not clear that disorder does not play a more 
crucial role. We will discuss the interplay of disorder and quantum critical behavior in both insulating and metallic systems in 
SectionESH 

There is of course no clear-cut way to separate systems whose behavior is governed by the proximity to a QCP and those in 
which disorder is the driving mechanism of NFL behavior A practical rule of thumb is to check whether the anomalous behavior 
exists in regions of the phase diagram distant from a clean quantu m ph ase transition, typically antiferromagnetism. This rule 
has been adopted, for instance, by Stewart in his review ( Stewart! 1200 ll) . Another indication of the importance of disorder in 
heavy fermion alloys is a large value of the residual resistivity = p{T ^ Q). Clean heavy fermion systems are characterized 
by large amounts of scattering o f the metallic carriers by the localized magnetic moments at high temperatures, the so-called 
Kondo scattering (Kondo, 1964). As the temperature is lowered, however, lattice translation invariance sets in and a precipitous 
drop of the resistivity, by orders of magnitude, occurs below the so-called coherence temperature Tcoh, which is typically a 
few tens of Kelvin. This delicate state is easily destroyed by disorder and large incoherent Kondo scattering is then able to 
survive to the lowest temperatures leading to values of po on the order of hundreds of pVl — cm. NFL behavior is observed in 
a number of heavy fermion alloys where no trace of coherence remains. In these cases, besides a large value of po, the leading 
low-temperature behavior is often linear in temperature with a negative c oefficient: p (T) ^ pn — AT, in sharp contrast to 
the Fermi liquid behavior of dilute Kondo impurities p (T) fa p^ — BT^ dNoziere si 11974 'Wilson'. 119751). Prominent exam- 
ples of systems where this linear in temperature behavior is observed are Yi -xPd.'^ lAndraka an d Tsveliklll99ltlOtF et alx 
1993: Seaman et ali fl99il) . UCus-^rPd^ (x = 1 and 1.5, se e Fig. |7j (Andraka and Stewart, 1993t IChau an d Maple. 1996 ; 
Weber ef qL.2001.) . UCu5_ ^Pt^ jx = 0.5, 0.75 and 1') dChau and Maple. .1996: .Stewart. .2001.) . Ui _-. Th^Pd9Alc cMaple ef 0/' . 
1995h . Ce j1.1Lan.3Cu2.2S i2 ( lAndrakalfT994l) . UCu4Ni dde la Torre et a/.LIl998h . URh^Geo (as grown) dSuUow ef a/.ll2000h . and 
U2PdSi2 (Li et ai, 1998). The last two systems and those of the form UCU4M, although nominally stoichiometric are actually 
crystallographically dis ordered, as evidenced by wide NMR lines (Bernal ef a/... .1995.) . pSR (.MacLaughlin et al... .1998i) and 
EXAFS measurements dBauer et al. l l2002HBooth et aail998il2002l) . 

Besides the anomalous transport, heavy fermion alloys also exhibit NFL thermodynamic properties. Singular behavior is 
observed in both the specific heat coefficient 7 (T) = C (T) /T and the magnetic susceptibility x (T) (see Fig.|8}. Ambiguities 
in the determination of the exact singularity are common since the available temperature range is often not very broad. Therefore, 
there are conflicting claims of logarithmic log (Tq/T)] or power law [~ divergences, a small value of A in the latter 

being very difficult to discern from the former dependence. Once again, we refer the reader to Stewart's review (Stewart, 200 ll 
for a summary of the published results. In that article, the author also presents many data previously published as being well 
described by a logarithmic dependence, replotted and fitted to a weak power law with a small A. These discrepancies highlight 
the difficulty of extracting the nature of the singular behavior from a limited range of temperatures. The only reliable way 
to resolve these ambiguities is through the determination of error bars for the exponents (as in classical critical phenomena). 
Unfortunately, this has not, to our knowledge, been attempted. 

Other experimental probes have also revealed NFL behavior in heavy fermion alloys. Optical conductivity studies have 
shown that the frequency dependence of the scattering rate at low temperatures is linear as opposed to the quadratic depen- 
dence predicted by Fe rmi liquid theory. This behavior has b een seen in UCu5_a;Pd2. (Degiorgi and Ott, 1996), Ux Yi_a;Pd3 
(iDegiorgi et al.l [19951) . and Ui_2;Th2:Pd2Al3 (iDegiorgi et all ri996) (see also the review (Degiorgi, 1999) for a more detailed 
discussion). Furthermore, neutron scattering data on UCu5_^Pdj; have revealed w/T scaling in the dynamical spin susceptibility 
over a broad tempera ture and frequency range without any significant wave vector dependence (other than the trivial Uranium 
atom form factor) ( Aron son et al.l 1995, 2001). Finally, we sh ould mention extensive NMR and pSR work on these systems 
jButtgen et aL .2000; .Liu et all 12000: MacLaughlin etal.l EoOL .1996> .2002a. .2004. .2002h. ^OOO. .20031 These studies, to- 
gether with the already cited EXAFS technique dBauer et al.l I2002t iBooth et al.l 1 19981 12002^ ■ are invaluable tools to access 
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Figure 7 Anomalous NFL resistivity of VCu^-xPdx- Data are from Ref. iAndraka and Steward 1 19931) . 
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Figure 8 Anomalous NFL magnetic susc eptibility of VCus-xPdx for different applied fields: H — 5 kOe (triangles) and H — 50 kOe 
(circles). Data are from Ref. iBemalef a/.LIl995.) . 



spatial fluctuations of local quantities, which are inevitably introduced by disorder. 



D. Disorder-driven NFL behavior in correlated systems: metallic glass phases 

1. Metallic glass phases 

Although usually not close to an ordered magnetic phase, disordered heavy fermion systems very often show spin-glass freez- 
ing. This is perhaps not too surprising given the interplay of disorder, local magnetic moments and RKKY-induced frustration 
in these systems. However, freezing temperatures are surprisingly low given the large concentrations of magnetic moments. A 
possible explanation is the Kondo compensation by the conduction electrons, which tends to favor a paramagnetic heavy-fermion 
state. Related to this is the proposal that the quantum critical point separating the heavy-fermion paramagnet and the spin glass 
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phase is at the origin of the NFL behavior of these alloys (see Section lTV.C.2> . Spin-glass behavior is usually detected by the 
difference between zero-field-cooled and field-cooled magnetic susceptibility curves. Direct measurements of the spin dynam- 
ics in /xSR experiments has a lso proven to be a usef ul tool. Examples of sys tems where spin-glass phase s have been detected 
are Y i.^^U^Pd.cj (x > 0.2 ) (iGaiewski et al.'. '1996'), UCus^^Pdj; (x > 1) ( fAndraka and Stewarl [l993; MacLaugMin et all 
1200 It IVollmer et all l2000l) URh^ Ge-^ CSullow ef a/.. 2000). UoP dSi-. llLi ef a/ 199&) . Ceo.isLao.ssCuiSia ( Andraka. 1994^ 
Uo.nr Thn a.^Ru^Si? llStewarll200ll Ui-^Y^Ab (0.3 <x< 0.7) llMavr ef adll997l) . and U2Cui7-xAl^ (x = 81 llPietri ef"^ 
Il997h . From these numerous examples, it becomes clear that spin-glass ordering is a common low-temperature fate of heavy 
fermion alloys. It is likely that a complete theoretical picture of the NFL behavior of these systems will have to encompass the 
spin-glass state. 

In many of these systems one can tune through a T = spin glass transition by varying parameters such as doping pressure, 
etc. In such cases one expects new physics associated with a spin-glass-paramagnet quantum critical point, but very few experi- 
ments exist where systematic studies of this regime have been reported in heavy fermion systems. On general grounds, one may 
expect such disorder-driven QCPs to have a more c omplic ated form then in the clean cases. In recent work, theoretical scenarios 
have been proposed ( Dobrosavlievic and Mirandal l2005t fVoita. 2003; Voita and Schmalian, 2004), suggesting that rare events 
and dissipation may "smear" such phase transition, possibly making it difficult even to detect the precise location of the critical 
point. More experimental work is urgently needed to settle these interesting issues, especially using the local probes (STM, 
NMR, /iSR) which are better suited to characterize such strongly inhomogeneous systems. 



2. Classiness in the charge sector 

An interesting alternative to spin glas s ordering is offered by t he pos sibility of glassy freezing of charge degrees of freedom. 
In many disordered electronic systems (iLee and Ram akrishnan. 1985), electron-electron interactions and disorder are equally 
important, and lead to a rich variety of behaviors which remain difficult do understand. Their competition often leads to the 
emergence of many metastable states and the resulting history-dependent glassy d ynamics of electrons. Theo retically, the 
possibility for glassy behavior in the charge sector was anticipated a long time ago fefros and S hklovskir. '1975') m situations 
where the electrons are strongly localized due to disorder In the opposite limit, for well delocalized electronic wave functions, 
one expects a single well defined ground state and absence of glassiness. The behavior in the intermediate region has proved 
more difficult to understand, and at present little is known as to the precise role and stability of the glassy phase close to the 
metal-insulator transition (MIT) (Mott, 1990). 

From the experimental p o int of view, glassy behav ior has often been observed in sufficiently insulating rn a terials 
(|pvadvahu^andP ollakl Il997t IVaknin et all Il998l l2000h . but more recent experiments jBogdanovich and Popovid l2002t 
Uaroszvhsk i et all l2002l) have provided striking and precise information on the regime closer to the metal-insulator transition. 
These experiments on low density electrons in silicon MOSFE T's have re vealed the existence of an intermediate metallic glass 
phase in low mobility (highly disordered) samples ( Bogdanovich and Popovic, 2002). Precise experimental studies of low tem- 
perature transport in this regime have revealed an unusual low temperature behavior of the resistivity of the form 

= p(o) + Ar3/2. 

This behavior is observed throughout the metallic glass phase, allowing one to systematically extrapolate the resistivity to T = 0, 
and thus characterize the approach to the metal-insulator transition. Such non-Fermi liquid behavior is consistent with some the- 
oretical predictions for metallic glasses ( Arrachea et al., 2004; Dalidovich and Dobrosavlievic, 2002; Sachdev and Read, 199ct), 
but the same exponent 3/2 is expe cted for glassiness both in the spin and the charge sector. To determine the origin of glassiness, 
the experiments were carried out jjaroszvnski et al 1 12004 in the presence of large magnetic fields parallel to the 2D channel 
(such that the field couples only to the electron spin), as shown in Fig |5] Since the Fermi energy is relatively small in this 
density regime (Tp ^ lOK), it is possible to completely spin polarize the electrons with accessible fields, and thus completely 
eliminate any spin fluctuations. Remarkably, all the signatures of the glassy behavior, as well as the T^/^ resistivity persisted in 
this regime, demonstrating that the anomalies originate from the charge degrees of freedom. 

These experimental results are certainly very compelling, as they indicate that glassy ordering in the charge sector may very 
well play a crucial role in most systems close to inhomogeneous insulating states. This possibility seems even more reasonable 
having in mind that screening effects can be expected to weaken as one approaches an insulator - thus strengthening the effects of 
the Coulomb interactions. More experimental work is necessary to test these ideas, especially using those techniques that directly 
couple to the charge degrees of freedom. In this respect, global and local compressibility measurements, as well as frequency- 
dependent dielectric constant experiments seem particularly promising. Such experiments have recently been reported on some 
transition-metal oxides, revealing glassy ordering, but similar work on other materials are called for in order to determine how 
general these phenomena may be. 
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Figure 9 Temperature dependence of the conductivity of two dimens ional electrons in silicon in the metallic glass phase, in the presence 
of a parallel magnetic field (B — 2 T). The data taken from Ref. ijaroszvnski et al\. |2004^ show results for several electron densities 
[ris (lO^'^cm^'^) = 11.9, 11.6, 11.3, 11.2, 11.0, 10.9, 10.7], as the metal -insulator transition is approached. 



IV. THEORETICAL APPROACHES 

A. Disordered Hubbard and Kondo lattice models 

1. Disordered-induced local moment formation 

The experiments of Paalanen et al. were very successfully described within a two-fluid model of itinerant carriers and localized 
magnetic moments (Paalanen et al., 1988; Ouirt and Marko, 1971), as explained in Section lTlLBI A first att empt to directly test 
this p henomenology within a well-defined model calculation was due to Milovanovic, Sachdev and Bhatt llMilovanovic et all 
Il989h . These authors investigated a disordered Hubbard model with both diagonal (site energies) and off-diagonal (hopping) 
disorder. The model parameters were chosen so as to faithfully describe the situation of Si:P. In contrast to the weak-disorder 
approach of the scaling theory, their work treats disorder exactly by numerical calculations while relying on the mean-field 
Hartree-Fock treatment of interactions, as was done for the single impurity Anderson model ( lAndersonlll96l . The latter is 
known to describe well the phenomenon of local moment formation, here taken to signify a broad temperature crossover with 
a Curie law magnetic susceptibility. However, it is not capable of correctly accounting for the low-temperature quenching of 
these moments by the conduction electrons, the Kondo effect. The results of Milovanovic et al. show that a finite fraction of 
the electrons do indeed exhibit a local moment Curie response even if the wave functions at the Fermi level are extended, i.e. 
the system is still metallic. The fraction of localized moments obtained is of order 10% at 25% above the critical density, in 
reasonable order of magnitude agreement with the experiments ( Paalanen et al., 1988). In other words, the results are compatible 
with the observation of local moment responses on the metallic side of the disorder-induced metal-insulator t ransition. A rather 
complete exploration of the phase diagram of this model with the same method has also been caiTied out jTusch and Loganl 
[1993), confirming this picture. 

Extremal statistics arguments were used by Bhatt an d Fisher to analyze loca l moment formation and also to determine the 
effects of residual interactions between these moments jBhatt and Fisheil|l992h . They considered a disordered Hubbard model 
very similar to the one studied by Milovanovic et al. In such a system, rare disorder fluctuations can give rise to sites weakly 
coupled to the rest of the lattice which, for sufficiently strong interactions, can form localized moments. In the absence of 
the Kond o effect or the RKKY int eractions between these moments, they would give rise to a Curie response as shown in 
reference jMilovanovic et a/1 Il989h . In the presence of the itinerant electron fluid in the metallic phase, however, the Kondo 
effect may quench these moments. The scale at which this quenching occurs is the Kondo temperature Tk ~ I?exp (— l/pj), 
where D is the conduction electron half band width, p is the density of states at the Fermi level, and J is the exchange coupling 
between the local moment and the conduction electron fluid. Spatial fluctuations lead to a distribution of couplings J and 
consequently to a distribution of Kondo temperatures. Thus, only those sites with Tk < T will contribute significantly to the 
thermodynamic properties. The overall behavior is still singular, though the corrections to the Curie law are very mild 

X {T) ^ ^ ^ i exp [-Aln'^ [In (To/T)]} , (40) 
where A and Tq are constants and d is the dimensionality. The effects of the RKKY interactions are stronger and lead, in the 
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absence of ordering, to a divergence that is slower than any power law 

x(r)~^~exp{Blni/'^(ro/T)}, (41) 

where again _B is a constant. The picture that emerges is that of a random singlet phase as will be explained below in Sec- 
tion llV.A^ A phase transition into a spin-glass phase, however, may intervene at low temperatures. 

This problem is further complicated by the fact that the Kondo temperature depends exponentially both on the exchange 
coupling J and the local conduction electron density of states p. The distribution of Tk calculated by Bhatt and Fisher did not 
take into account the fluctuations of the latter. We emphasize that p should be generalized, in a disordered system, to the local 
density of states p (R) at the position R of the magnetic moment. This is not a self-averaging quantity and its fluctuations are 
expected to grow exponentially as the dis order-induced metal-insulat or transition is approached. These effects were analyzed by 
Dobrosavljevic, Kirkpatrick and Kotliar (iDobrosavhevic et aZ.[ll99 j) . These authors emphasized the universal log-normal form 
of the distribution of local density of states, which in turn leads to a universal Tk distribution. As usual for a broad distribution, 
the magnetic susceptibility can be written as x {T) ^ n (T) /T, where n (T) ~ T^C^) is the temperature dependent number 
of unquenched spins at temperature T. Although the exponent a (T) ^ as T — » 0, it does so in an extremely slow fashion 
and in practice, for reasonable values of the parameters and for a strongly disordered but still metallic system, n (T) ~ 50% for 
T - 10-4 - IK. 



2. Random singlet phases 

Local moment formation is generically a crossover phenomenon which occurs at intermediate temperatures whose signature 
is a Curie-Weiss magnetic susceptibility. Possible low-temperature fates of these moments include some form of magnetic 
ordering, such as antiferromagnetism and ferromagnetism or spin-glass freezing. A novel kind of low temperature fixed point of 
disordered localized magnetic moments, ho wever, was proposed from s tudies of doped semiconductors close to the localization 
transition, as was discussed in Section lTtl.BI I Bha tt and Leelll98ll [l982^. In these systems, the magnetic susceptibility shows no 
signs of any type of ordering but diverges with a non- trivial power law CPaalanen and Bhatt,. 1 99 k.Sarachik.. 1995.) 

X{T) - (42) 

7(T) = ^ ^ r"-i. (43) 

The theory of Bhatt and Lee ("Bhatt and Lee, '198 1', '1982') was based upon a generalization to three dimensions of a method 
introduced by M a, Dasgupta and Hu for one-dimensional disordered antiferromagnetic spin chains (.Dasgupta and k. Ma, 198(i 
iMa^f fl/.Lfl979l) . The Hamiltonian of such a random spin system is written as 

iJ^^Jy S, -Sj, (44) 

where > is a random variable distributed according to some given distribution function Pq (Jij). The method consists of 
looking for the strongest coupling of a given realization of the random lattice, say fl. At energy scales much smaller than ft, 
there is only a small probabiUty for this pair of spins to be in its excited triplet state. We can then assume the pair is locked in 
its ground singlet state and is magnetically inert at the scale considered. The pair is then effectively "removed" from the system. 
The remaining spins which had interactions with these locked spins develop additional interactions induced by the removed pair 
If spins i and j form the removed pair and spins k and I interact with i and j, respectively, the new interaction can be obtained 
in second order of perturbation theory 

AA, = <45) 

If initially k and I do not interact they will do so after the decimation step. Their new coupling will be smaller than any of the 
removed ones. Therefore, as more and more spins are decimated, the largest couplings are progressively removed while new 
smaller ones are generated, thus changing the shape of the "bare" distribution into an effective "renormalized" one at lower 
energy scales. Formally, the renormalization group flow starts at P [J, Hq] — Pq ( J), where flo is the initial largest value of 
the distribution. After many decimations, fto will have been decreased to ft and the distribution will be P [J, ft]. Bhatt and 
Lee implemented numerically the decimation procedure we have just outlined using realistic initial distributions of couplings 
appropriate for the insulating p hase of doped semiconducto rs. They showed that these initial distributions flowed towards very 
broad distributions at small J (iBhatt and Leei Il98l[ Il98l . leading to diverging thermodynamic responses of the form given 
in Eq. ( I43> , in remarkable agreement with the experiments. The low energy state of the system showed no tendency towards 
ordering but was characterized by the successive formation of singlet pairs between widely separated spins at each decreasing 
energy scale. This novel disordered magnetic state was dubbed a "random singlet phase". 
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3. Phenomenological Kondo disorder model 

Although the importance of the distribution of Kondo temperatures was first proposed with doped semiconductors in mind, 
the context of disordered heavy fermion materials, where the Kondo effect is of primary importance, seemed like a natural 
arena for its applications. Indeed, it was put to good use in the attempts to understand both the temperature dependence of 
the Cu NMR line-w idths and the thermodynamic properties of the Kondo alloy UCu5_j;Pda; {x ~ 1, 1.5) by Bernal et al. 
( iBernalef adll995h . As was mentioned in Section llll.CI this compound is known to exhibit non-Fermi liquid behavior in many 
of its properties. Bernal et al. analyzed the broad Cu NMR lines in the following way. The Knight shift at a particular Cu site 
positioned at R, K (R), is proportional to the local spin susceptibility 

i^(R)=a(R)x(R), (46) 

where a (R) is the hyperfine coupling. If a (R) is assumed to have little variation in the sample, the spatial fluctuations of the 
Knight shift can then be ascribed to the spatial variations of the local susceptibility 5K = aSx- If we use the following fairly 
accurate form for the susceptibility of a single Kondo impurity 

x{T)- ^2 rr ' (47) 

where as usual Tk ~ I?exp [—1/ {pJ)], we can use the distribution of Kondo temperatures P (Tk) to find the variance of the 
susceptibility and hence the line-width. Furthermore, the bulk susceptibility and the specific heat can be obtained by performing 
the same kind of average procedure over single-impurity results. Bernal et al. first fitted the susceptibility data by adjusting 
the mean and the variance of a Gaussian distribution of the bare quantity pJ. Then, a fairly reasonable agreement with the 
NMR line-width and specific heat data were obtained without further adjustments, explaining in particular the strong anomalous 
temperature dependence of the line-widths. The distribution of pJ used was not too broad. However, due to the exponential 
formofT/f ~ D exp [— 1/ (pJ)], even a fairly narrow distribution of p J can lead to a very broad P (Tk). In the case of Bernal 
et al., P (Tk) — > const, as Ti<- —> (see Fig.[T0j, so that 



P{Tk) 



T + aTK 







dTK ^ h const. 

T + oTk 



- In (I) , (48) 

in good agreement with the experimental findings (see the solid lines of Fig|8}. The picture is again one where a few spins whose 
Tk < T (shaded area in Fig.llQ> dominate the thermodynamic response. 



4. Dynamical mean field theory (DIVIFT) 

The phenomenological model of Bernal et al. made the seemingly unjustified assumption of a collection of independent single 
Kondo impurities, even though the alloys UCu5_j.Pd^ contained a concentrated fee lattice of U ions. This assumption was put 
on a firmer theoretical basis by Miranda, Dobrosavljevic, and Kotliar (Miranda ef a/., 1996, 1997a b). These authors used a 
dynamical mean field theory (DMFT) approach to the disordered Kondo (or Anderson) lattice problem. The DMFT method has 
become a very useful tool in the study of strongly correlated materials. It is the natural analogue of the more familiar Weiss 
mean field theory of magnetic systems, whic h is here generalized to describe quantum ( especially fermionic) particles. Fairly 
complete reviews of the subject are available (iGeorges et a/.LfT996tlPruschke et a/.lll995l) . and we will content ourselves with a 
brief description, emphasizing the physical aspects. 

For definiteness, let us focus on a disordered Anderson lattice Hamiltonian in usual notation 

Hand ^ '^Y^ {^a^ja + H- c-) + Yl ^jS^S^ 

{ij)a jcr 

io- j 
+ E(^4<.S-+H.c.), (49) 
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Figure 10 Distribution of Kond o temperatures obt ained within the Kondo disorder model from an analysis of the disordered heavy fermion 
compounds UCus-^Pdj, (Ref. iBemal et a/.LIlQQ.^ ). Spins with Tk < T (shaded area) dominate the thermodynamic response. 



where the conduction, /-electron, and hybridization energies e^, Efj, and Vj are, in principle, random variables distributed 
according to (e^), Pe^ (Efj), and Py (Vj), respectively. As in the Weiss theory, one starts by singling out a particular lattice 
site, say j, and writing out its effective action. A simplification is made here, by neglecting all non-quadratic terms in the local 
fermionic operators, except for the [/-term. One gets 
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The site j "talks" to the rest of the lattice only through the bath (or "cavity") function Acj (r) in Eq. ( I51> . For simplicity, we 
particularize the formulation to the case of a Bethe lattice, in which the bath function is given by 



A, (r)=t^G^-(r). 



(54) 



Here, the average local Green's function is obtained by first calculating the local conduction electron Green's functions governed 
by the actions ( 15 0> 



V w4(o) 



eff ' 



(55) 



and then averaging over all sites with the distributions of bare parameters P^, Pej, and Py. Thus, all local coiTelations are 
accounted for faithfully, while inter-site ones enter only through the bath function. As in the Weiss theory, the whole procedure 
can be shown to be exact in the limit of large coordination ("infinite dimensions") if appropriate scaling of the hopping is 
performed {t ^ t/^, with t held constant as z ^ oo). Unlike the Weiss theory, however, the "order parameter" here, 
namely, Ac (t) is a. function not just a number (hence the name "dynamical"). Correlation effects are incorporated in the step 
where the conduction electron Green's function (I55> is calculated. This is equivalent to solving a single-impurity Anderson 
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model jGeorges and Kotl iar. '1992') and since the bare parameters are random, one has to solve a whole ensemb le of these. Th e 
treatment of disorder within DMFT is equivalent to the well-known coherent potential approximation (CPA) ( Elliott et a/lll974l) . 

The connection to the Kondo disorder model can now be made more apparent. The ensemble of single impurity problems 
described by the actions (15 0> is the equivalent of the collection of independent single Kondo impurities of the phenomenological 
approach. However, the single impurity problems are not really independent as each one "sees" the same bath function ( I54> . 
which in turn contains information from all the other sites. The fully self-consistent calculation of the bath function, however, 
shows that its precise form does not contain essential features being as it is an average over many different sites. Each Anderson 
impurity has its local Kondo temperature Tkj, which in the Kondo limit can be written as (taking U oo) 

Tk, = D exp [- \Efj I / (2pF/)] . (56) 

The connection to the Kondo model is obtained from Jj — 'iVf I \Efj\. A distribution of Kondo temperatures then follows. 
The full solution is able to produce, with a judicious choice of bare parameters, a P (Tk) very similar to the experimentally 
determined one for UCu5_a;Pda; based on the simple phenomenological Kondo disorder model (Section riV.A.3> . 

Within the DMFT framework, one can proceed to the calculation of various physical properties. Bulk thermodynamic re- 
sponses in the particular case of the Anders on/Kondo lattice, can be very accurately obtained through an ensemble average of 
the individual contributions from each site (iMiranda et all Il996h . thus justifying the procedu re adopted by Bernal et al. . The 
reason for this is the dominance of the /-site contributions over the conduction electron part (iMiranda et al.l fT99^. Transport 
properties, however, cannot be so easily calculated and require a special consideration of the coherence of the motion through 
the lattice (Tesanovic, 1986). This is particularly striking in the clean limit, in which it is well known that the /-ion contributions 
very accurately add up to give the thermodynamic properties but the onset of heavy fermion coherence in the resistivity is totally 
absent from a single impurity description. 

It is precisely with regard to transport properties that the DMFT approach is able to go beyond the phenomenological Kondo 
disorder model. Assuming the same Tk distribution obtained experimentally by Bernal et al ., one gets a resistivity wh ich 
has a large residual value and decreases linearly with te mperature with a nega t ive sl ope ( Mir anda et all 1 19961 Il997ald) . in 
complete agreement with the experiments on UCu5_^Pdj; ( Andraka and Stewart, 1993). The large residual value is due to the 
total destruction of coherence by disorder The anomalous NFL linear dependence comes from the gradual unquenching of local 
moments as the temperature is raised. 

Other physical properties were calculated within DMFT, with good agreement with experiments: the dynamical spin sus- 
ceptibility ( Aronso n ef g/., 1 995; Miranda ef a/., 1996)^ the jnagneto-resistance (Chattopadhvav ef a/., 1998), and the optical 
conductivity dChattopadhvav and Jarrelllfl997tlDegiorgi and Otttll996l) . In all cases, NFL behavior was tied to the finite weight 
of P{Tk) as Tk 0. 

One of the perceived deficiencies of this theory is its extreme sensitivity to the choice of bare distributions. This is due to 
the exponential dependence of the Kondo temperature on model parameters, which is hardly affected by the self-consistency 
introduced by DMFT. In the particular case of UCu5_2;Pda;, the finite intercept P {Tk — 0) is the most important feature of the 
distribution and only fine tuning of the bare distributions can provide the c orrect value. Besides, the thermodynam ic properties 
of many other compounds are more accurately described by power laws ( de Andrade et al.l Il998t Istewarti EoOll) . These can 
be accommodated within the theory through a power-law distribution of Kondo temperatures P (Tk) ~ However, once 

again fine-tuning is necessary if one wants to obtain such a distribution within DMFT. 

5. Statistical DMFT: localization effects and results of numerical calculations 

The problem of fine tuning within DMFT was remedied in a more complete approach. An important assumption of DMFT is 
the averaging procedure contained in Eq. (|53- It means that each site feels an average hybridization with its neighbors, which 
is why the procedure becomes increasingly more accurate the larger the coordination. Fluctuations of this quantity, usually 
associated with And erson localization effects, can, however, be incorporated in a so-called Statistical Dynamical Mean Field 
Theory (statDMFT) ("Dobros avhevic and G. Kotliar. '1997'; Dobrosa vhevic and Kotliail '1998'). A discussion as applied to the 
Anderson lattice can be found in reference ( Aguiar ef a/., 2003). The assumption of retaining only on-site correlations, as in 
DMFT, is maintained in statDMFT. Formally, this means that one still has to solve an ensemble of single impurity problems 
defined by the actions (15 0> . The only difference comes from the prescription for the bath function, which is now site-dependent: 
Ac (r) Ac {j, t) . If the lattice coordination is z, we can write 

z 

Ac(j,r)=i2^G«(/,T), (57) 
i=\ 

where G^P {I, r) is the local Green's function of nearest-neighbor site I, after the removal of site j. Once again, we are partic- 
ularizing to the Bethe lattice, which simplifies the discussion considerably. The Green's function with a nearest-neighbor site 
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removed is defined as in Eq. i55\ but has to be calculated with a different action. This action has the same form as ( I50t . but now 
the bath function Ac (r) Ac^ {I, r), where 

z-l 
m—l 

Note that in (I58> the summation has only z — 1 terms, as opposed to the z terms of i57\ . since site j has been removed. The 
self-consistency loop is closed by noting that the Green's functions on the right-hand side of Eq. ( I58> are the same (i. e., 
obey the same distributions) as the ones on the right-hand side of Eq. (I57t . We emphasize that the self-consistency problem 
is now no longer a simple algebraic equation but rather a set of stochastic equations in the functions Gc (j, t) and ci"'' (I, r), 
whose distributions have to be determined. If the inte ractions are turned off, the treatment of disorder we have described is 
equivalent to the self-consistent theory of localization llAbou-Chacra et a7l Il973h . The latter is known to exhibit a disorder- 
induced Anderson metal-insulator transition for coordination z > 3. The formulation above can be generalized to any lattice 
fDobrosavIievic and Kotliar, 1998). However, the full set of stochastic equations has to be solved numerically, in which case 
it is advantageous to work in a Bethe lattice. In a realistic lattice, this approach offers the advantage of treating the disorder 
exactly, albeit numerically, while incorpora ting the local effects of correlations. In this sense, it is a natural generalization of 
the work of Milovanovic et al. ( Milovanovic et al. 1 I1989I) . to which it reduces if the single-impurity problems are treated within 
Anderson's mean field theory (Anderson, 1961'). More sophisticated treatments of the single-impurity problems are, however, 
possible, especially in the low temperature region. 

The ability to access full distribution functions of physical quantities is a especially appealing feature of statDMFT. In partic- 
ular, the local density of states 

Pe(i,c^) = -Im[GeO>-i<5)] (59) 

TT 

is intimately connected with the transport properties. It has the natural interpretation of a escape rate from site j and is expected 
to be finite if states at energy lo are extended and to vanish if they are localized, as first pointed out by Anderson in his ground- 
breaking work ( Anderson. 1958.) . Furthermore, the distribution of can become very broad with increasing disorder 
Although the average density of states pc (w) remains finite for any disorder str ength, its typical value p^J'P (uj) vanishes at 
the metal-insulator transition and can serve as an order parameter for localization dAndersonl 1 1 95 8h . A convenient measure of 

pI^p (uj) is the geometric average cxp |ln [p 

Furthermore, other local quantities of interest can also be computed, such as the distribution of Kondo temperatures P (Tk) 
already familiar from DMFT. Recalling that Tk ~ Dexp [—1/ (pJ)], we notice that the density of states p "seen" by any 
particular site j is a fluctuating quantity in statDMFT and is thus a new source of Tk fluctuations. For the case of the /-electrons 
of the disordered Anderson lattice, for example 

T ITT- ^ ■ (60) 

In particular, even if J = 21/^/ \Ef \ is not random (e. g., if only conduction electron diagonal disorder is present), there will be 
a distribution of Kondo temperatures. Moreover, in DMFT, a discrete distribution of J leads necessarily to a discrete distribution 
of Tk, as p is fixed. This is no longer true in statDMFT, however. Even if J (or, equivalently, V or Ef) is discrete, there will 
be a continuous distribution of p and consequently of Tk- This is because pj at a particular site is influenced by fluctuations of 
other quantities at sites which may be very distant from j, due to the extended nature of the conduction electron wave function 
in the metallic phase. 



Pj = — Im 



6. Electronic Griffiths pliase in disordered Kondo lattice models of dirty heavy-fermion materials 

The application of statDMFT to the disordered Anderson lattice problem was described in references ("Aguiar et al l l2003t 
[Miranda and Dobrosavlievic, 1999, 2001a b). The single-impurity problems were solved at T = within the slave boson mean 
field fi eld theor y (Barnes, 1976; Coleman, 1984, 1987; Read andNewns, 1983) or for T ^ {) with second-order perturbation 
theory dKaiueter and KotliaL.1996;.Mever and Noltina. |2000a b). Several important features deserve mention. The distribution 
of Kondo temperatures is genericaUy log-normal for weak disorder but slowly evolves towards a power-law at small Tk, 

P[Tk)^T^-\ (61) 

where the exponent a is a continuously varying function of the disorder strength W , as shown in Fig. ^2 This generic form 
is fairly insensitive to the particular form of the bare distribution. This universality reflects the mixing of many single-impurity 
problems which are connected by the extended wave function of the conduction electrons within a correlation volume. 



26 



1.0 



0.0 



?5 -1.0 
O -2.0 



-3.0 



-4.0 




W/t=0.5 

W/t=1.0 

W/t=1 .5 

W/t=2.0 

W/t=2.5 



-12.5 



-7.5 -2.5 

x=log(TJ 



Figure 1 1 Distribution of Kondo temperatures within statDFMT, showing P (Tk 



strength of disorder. The onset of NFL behavior occurs at a < 1. From Ref. iMiranda and Dobrosavlievid.l2001bl) . 



with g (ty) con tinuously varying with the 



From the power-law distribution of energy scales, the usual phenomenology of a Griffiths phase follows. Although our dis- 
cussion in this Section is self-contained, a complete description of a generic quantum Griffiths phase is deferred to Section UV BI 
From the distribution ibH . the susceptibility and the specific heat coefficient, for example, are given by 



xiT) 



C{T) 
T 



dTK 



P{Tk) 
T + aTx 



(62) 



For sufficient disorder, a < 1 and NFL behavior is observed. The marginal case a — \ corresponds to a log-divergence. Note 
that the point a = \ does not signal any phase transition. Other physical quantities such as the non-linear susceptibility x^'^'' {T) 
start to diverge at a larger value of a (corresponding to smaller disorder strength), since x'^' (0) ~ 1 /T^. Several other physical 
quantities can be obtained as continuous functions of the exponen t a. Such behavior appears consistent w ith many observations 
in disordered heavy fermion materials which show NFL behavior dde Andrade etal l ll998l:IStewarll200ll) (Section HlTci . 

Griffiths phases are common in the proximity of magnetic phase transitions, as is explained in more detail in Section lTV.BI 
Here, however, the Griffiths phase occurs in the absence of any form of magnetic ordering. Rather, it is associated with the 
proximity to the disorder-induced metal-insulator transition, which occurs for sufficiently large disorder (Wmit ~ 12t for 
conduction electron diagonal disorder only (Miranda and Dobrosavlievic, 2001b)). This is somewhat reminiscent of the local 
moment phase close to the metal-insulator transition in doped semiconductors, but is different in two important respects: (i) 
in a disordered Kondo/Anderson lattice, local moments are assumed stable even in the absence of disorder and (ii) the NFL 
divergences begin to show up at fairly weak randomness and persist for a wide range of disorder strength. Many other features 
of the statDMFT solution of the disordered Anderson lattice, particularly concerning transport properties, can be found in 
references ( Agtliar et al., 2003; Miranda and Dobrosavlievic, 1999, 2001a b). 

We should stress that the anomalous power laws obtained in this approach are a direct consequence of the power law in the 
distribution of Kondo temperatures. Interestingly, all that is needed is a power-law distribution of energy scales for the spin 
fluctuators, not necessarily of Kondo origin or related to localization effects. Given this form of distribution the associated 
Griffiths divergences follow immediately (see a discussion of generic properties of Griffiths phases in Section IIV.B.2> . As a 
result, widely different microscopic mechanisms can give rise to the same macroscopic behavior. This will be exemplified later 
when we discuss magnetic Griffiths phases (Section IIV.B> . which are usually tied to a magnetic phase transition. Whereas 
this highlights an interesting universality of this phenomenology, it makes it hard to distinguish between different microscopic 
mechanisms based solely upon macroscopic measurements. 

Some insight into the origin of the power law distribution of Tk was offered in reference (iTanaskovic et a/.Ll2004) . In that 
work, it was shown how a power-law distribution of Kondo temperatures can be easily obtained within a DMFT treatment of the 
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Figu re 12 Comparison bet ween the a (W) exponent determined within statDMFT and the one obtained in the effective model of 
Ref. jTanaskovic et a/.Ll2004) . 

disordered Anderson lattice, if a Gaussian distribution of conduction electron on-site energies 



is used, with fixed values of Ef and V. In this case, although the bath function is fixed at its average value (Section [lV.A.4t . 
different values of ej generate different Kondo temperatures (cf. Eq. i50\ ) according to 



where ^4 is a constant and T]^ is the Kondo temperature for ej — 0. It is then straightforward to show that, up to log corrections. 



where = 2AW'^. This type of argument, in which an energy scale depends exponentially on a random variable, which 
in turn occurs with an exponential probability, leading to a power law distribution of the energy scale is very common in, and 
perhaps generic to Griffiths phases (see Section llV.B.2> . Now, although the specific form of the distribution i63\ is required 
within DMFT for the power law to emerge, in statDMFT fluctuations of the bath function are generically gaussian, at least at 
weak disorder By noting that the bath function enters additively to tj in the single-site action (cf. Eq. ( 15 1» . it is clear that 
the appropriate statistics is naturally generated within this approach. Careful comparison between the two approaches confirms 
that this is indeed the case (see Fig.ll2>. thus clarifying the generic nature of the mechanism behind the Griffiths phase in these 
systems. 

7. Electronic Griffiths phase near the IVIott-Anderson metal-insulator transition 

The StatDMFT approach can be usefully employed in the analysis of the disordered Hubbard model, which serves as a proto- 
typical model for both disorder-induced (Anderson) as well as interaction-induced (Mott-Hubbard) metal-insulator transitions. 
The Hamiltonian in this case is 
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Figu re 13 Distribution of Kondo tem peratures (here named Z) within statDMFT for the disordered Hubbard model. From 
Ref. iPobrosavlievic and G.KotMa3.ll997l) . 

The StatDMFT equations in this case are very similar to the Anderson lattice and we write them down here for completion. The 
auxiliary single-site actions read (.Dobrosavlievic and G.Kotliar. 1997: Dobrosavlievic and Kotliar„.1998ii) 

S^s'' ij) f dr'cl (r) [S (r - r') 

^ Jo Jq 

X (dr + e, - m) + Ac (j, r - r')] c^, (r') 

+ U [ drc]^ (t) c^^ (t) c]^ (r) c^.^ (r) , (67) 
Jo 

where the bath function is once again given by Eq. ( I57> . The Green's function with a nearest-neighbor site removed is obtained 
from an action of the form ( I67> . with a bath function given by Eq. (I58> . which closes the self-consistency loop (again, in a Bethe 
lattice). 

The analysis of the distributions of two local quantities which come out of the single-site actions \bl\ provide especially useful 
insights into the physics of the disordered Hubbard model. One is the local density of states at the Fermi level p (j, 0), already 
defined in Eq. ( I59> . The other is the local Kondo temperature of each single-site action Tkj. The product of these two is the 
quasi-particle weight Zi = p (j, 0) Tkj- While the Kondo temperature governs the thermodynamic response of the system, the 
local density of states is related to the transport properties, as explained in Section lTV.A.SI 

The results of statDMFT as applied to the disordered Hubbard model show that, similarl y to the disordered Anderson lattice. 
a Griffiths phase emerges, w ith NFL behavior, for intermediate values of disorder strength jDobrosavlievic and G.Kotlia3.ll997t 
iDobrosavlievic and Kotliatj 0*998) . This is signaled by a power-law distribution of Kondo temperatures P (Tk) ~ where 
a varies continuously with disorder. The marginal value of a = 1, at which the magnetic susceptibility diverges logarithmically 
with decreasing temperature, occurs atW 7t, as shown in Fig.[0] This is reminiscent of the behavior of doped semiconduc- 
tors. In contrast to the Anderson lattice though, the NFL behavior requires much stronger disorder. This can be rationalized by 
noting that in contrast to the latter, where local moments are stable in the clean limit, here they require quite a large amount of 
randomness to appear. 

Moreover, for low conduction electron filling (n w 0.3), the typical local densi ty of states goes to zero at a large 
value of disorder strength Wm it ~ Hi, signaling a metal-insulator transition (iDobrosavhevic and G.Kotlia3 Il997t 
iDobrosavhevic and Kotliatj Il998h . As in the non-interacting case, the average value remains non-critical, actually diverging 
at Wmit- The full distribution of p {j, 0) is very close to a log-normal form, with a width which increases as disorder increases. 
Closer to half-filling {n k, 0.7), however, the behavior is quite different. At these larger densit ies the typical value of p (7, Q ) 
remains finite for very large values of disorder, suggesting that Wc is pushed up considerably (.Dobrosavlievic and KotliaJll998h . 
A similar "scree ning" of the disorder po tential by strong interactions has been also found in a DMFT study of the disordered 
Hubbard model (iTanaskovic et al\\200'^ . 

8. Incoherent metallic phase and the anomalous resistivity drop near the 2D MIT 

The puzzle of the two-dimensional metal-insulator transition observed in several materials ( lAbrahams et al\. 1200 ll) remains 
a largely unsolved problem. One of the most intriguing features of these systems is the large resistivity drop observed in a 
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temperature range comparable to the Fermi energy. This suggests the importance of inelastic scattering processes, a feature 
usually absent in Fermi liquid based approaches, or at most included only in a perturbat ive fashion. An initial attempt at the 
incorporation of such inelastic effects was made within a DMFT framework ( Agu iar et al. 112004). A half-filled disordered Hub- 
bard model was considered in that work. The associated single-impurity problems were solved within second-order perturbation 
theory cKaiueter and Kotliar . 1996) and checked against Quantum Monte Carlo calculations. The results were also compared 
with a Hartree-Fock impurity solver, which does not include inelastic scattering. 

The most important result of that work is the demonstration of the large contribution of inelastic processes to transport 
properties in a region of parameter space where kinetic energy, interaction and disorder strength are of comparable size. This 
is especially striking when one compares the results of DMFT with Hartree-Fock. Whereas in the former, a drop of order 
10 in the total scattering rate is obtained, only a very weak drop is observed in the latter for the same values of disorder and 
interaction. This large difference is clarified when one separates the inelastic from the elastic parts. The Hartree-Fock scattering 
rate, which is completely elastic, is comparable to the elastic contribution in DMFT. However, for comparable values of disorder, 
Fermi energy and interaction strengths, the contribution of inelastic scattering is seen to completely overwhelm the weak elastic 
processes, in effect determining the size and temperature dependence of the total rate. 

Furthermore, keeping the interaction and Fermi energy of equal size and decreasing randomness, one is able to enhance even 
more the size of the drop in the scattering rate. The more gradual drop for larger disorder is attributed to a wide distribution 
of coherence scales. These scales set the boundary between coherent transport at low temperatures and incoherent inelastic- 
process-dominated transport at higher temperatures. Once a distribution of these scales is generated, different spatial regions 
of the system become incoherent at different temperatures as the temperatures is raised, leading to a gradual rise in the overall 
scattering rate. This rise is approximately linear in the range where the various coherence sca les are equally frequent. An 
approximately linear rise of the scattering rate is observed in Si MOSFET's ( lAbrahams et a/.Ll2001 ). Finally, the coherence scale 
is inversely proportional to the local effective mass, which in turn is tied to the thermodynamic response of the system. Thus, 
when the disorder is comparable to the interaction strength, a sufficiently broad distribution of coherence scales is generated such 
that very small scales have finite weight, leading to an enhanced overall effective mass and thermodynamic properties. Results 
similar to those DMFT findings we re also obtained from exact numerical studies of finite size lattices ( Denteneer and ScalettaJ 
l2003tlDenteneer et a/.LIl999[l200ll) . but more work is needed to gain a more definitive understanding of this regime. 

B. Magnetic Griffiths phases 

1 . Quantum Griffiths phases in insulating magnets with disorder: random singlet formation and the infinite-randomness fixed 
point 

We have seen in Section IIV.A.2I how the Bhatt and Lee numerical study of random antiferromagnetic spin systems 
(^hatt and Lee, 1981, 1982) led to the concept of a random singlet phase. Further insight into this kind of behavior came 
from a rather complete analytical treatment of the nearest-neighbor random antiferromagnetic chain (Eq. \AA\ ) by D. S. Fisher 
(lFisheilll994l) using the renormalization group procedure invented by Ma, Dasgupta and Hu and described in Section ITV. A. 21 
Fisher showed that an essentially exact solution of the problem could be obtained asymptotically at low energies. He showed 
that the relevant fixed point distribution of couplings is infinitely broad and the system slowly approaches it as it is probed at 
lower and lower energies. The system is governed by an "infinite randomness fixed point" (IRFP) which reflects the true phys- 
ical significance of the random singlet phase. Note that, if the distribution is very broad, the side couplings of the strongest 
bond il, Jik and Jji in Eq. ( 145 > . are almost certainly much smaller than VL, rendering the perturbation theory result essentially 
exact. It is in this sense that the Fisher solution is asymptotically exact. Initial decimations will be inaccurate but if the system 
flows towards broader and broader distributi ons, then the results wil l become increasingly more accurate and asymptotically 
exact. Interestingly, results at weak disorder llDotv and Fished 1 1992 1) showed that the clean system is unstable with respect to 
infinitesimal disorder and numerical results seem to confirm that practically any initial randomness is in the basin of attraction 
of the IRFP found by Fisher We will highlight the some of the most important physical properties of this type of phase. This 
is a vast subject and we r efer the reader to some recent reviews of experimental and theoretical results of disordered spin chains 
and ladders ( Continentin o et gP. 12004: Igloi and Monthus, 2005). 

When the system is governed by the IRFP it exhibits many unusual properties. In particular, the relation between energy (il) 
and length (L) scales is "activated" 

17-17oexp(~L'^), (68) 

where ijj = 1/2 in the Heisenberg chain case. This should be contrasted with the more common dynamical scaling observed in 
quantum critical systems Q, ~ L^^', where z is the dynamical exponent (Hertz., 1976) (see Section ITlI. A. H . Formally, we can 
say that the "activated dynamical scaling" of Eq. ( I68> corresponds to a divergent dynamical exponent z ^ oo. This has many 
important consequences. The magnetic susceptibility of the system at temperature T can be obtained by scaling down to f2 = T. 
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It can then be very accurately calculated through 



(69) 



where it is assumed that the spins that were decimated at larger scales form magnetically inert singlets while the remaining ones, 
whose density is n (il = T) are almost free and contribute a Curie-like term. This estimate is better, the broader the distribution 
of couplings. If we now write 

{n = T)^L{n^T)^ [In iflo/T)]''/''' , (70) 
where use was made of Eq. ( I68> . we can write for the random Heisenberg chain 

^^^^^rini/'/'(^) ^T\n'{^)- ^^^^ 
Similarly, the entropy density s (T) ~ n (fi = T) In 2, from which follows 

C (T) IT ~ ^ ~ ^ (72) 

Furthermore, the spin-spin correlation functions are very broadly distributed. Fisher showed that while the average correlation 
function decays as a power law (Cij = (S^ • Sj)) 

- ^^^^ 

\^-J\ 

the typical one, given by the geometric average, decays as a stretched exponential 

exp (h[\C~\) ^ exp (- K - J l'^) ^ exp (- - ]{) . (74) 

This large difference between average and typical correlations is a result of the fact that a typical pair of spins separated by a 
distance L ^ 1* — j | does not form a random singlet and is therefore very weakly correlated. Those rare pairs that do form a 
random singlet, however, are strongly correlated and dominate the average. 

Further extensions of these studies deserve mention. Random chains with highe r spins have been investigated. Integer spin 
chains are especially interesting because the clean system is gapped (Haldane gap ('HaldaneVl983a'b')). It has been shown that 
the S — 1 chai n (and pres umabl y othe r integer S ones) are stable with respect to weak eno ugh disorder, r etaining a (pseudo- 
)gap (Boechat ejM Il99fit lHidaill999t |H^an and Yang, 1997; Monthuj.eI^ 1 19971 Il998t ISaguia a/.L l2003atl Yang et al} 
ll996l) For larger disorder, the pseudogap is destroyed and a Griffiths phase is realized, with a diverging power-law susceptibility 
with a varying non-universal exponent. At sufficiently large disorder, a random singlet phase is obtained governed by an IRFP. 
Similar results were also obtained in other sy stems for which the c lean anal ogue is gapped: random chains in which the gap 
is induced by a dimerization of the couplin gs (iHeneliu s and Girvinl 11998k: .Hvma n et a/.l [19961 lYang et a/.LfT99 6) and the two- 
leg ladder (Hovos and Miranda, 2004; Meli n et aZT l2002t lYusuf and Yangl l2002l) . A Griffiths phase was also identified in a 
disordered Kondo necklace model (Doniach, 1977), in which the spin sector of the one-dimensional Kondo lattice is mimicked 
by a double spin chain system (Rappoport et ai, 2003) . 

It has been proposed that the = 3/2 antiferrorn agnetic random chain realizes two kinds of random singlet phases separated 
by a quantum phase transition (iRefael ef a/.ll2002t) . At weaker randomness, the excitations of the random singlet phase have 
5=1/2 whereas in the stronger disorder phase they have S = 3/2. The quantum critical point between these two phases 
can be vie wed as a multicritical p oint in a larger parameter space which includes the effects of an added dimerization of the 
couplings dPamle and Husel|2002ft . The physics of this and other higher spin multicritical points has been shown to be governed 
by infinite randomness fixed points with different exponents: for example, i/; — 1/4 for the S = 3/2 chain. T he phase diagram of 
the S = 3/2 chain proposed in Ref. dRefael et a/.Ll2002h has been recently disputed, however dCarlon et a/.t(2004jilSaguia et all 
I2003bh . 

Another important universality class of disordered spin chains is found in systems with random ferromagnetic and antiferro- 
magnetic interactions. In this case, it is possible to use a generalization of the de cimation procedure describe d above for any 
spin size a nd both signs of cou plings (Frischmuth and Si^rist. 1997: Frischmuth et ai . 1999: Hida. 1997: Hikihara et a/.[ll999t 
IWesterber g et al l fT99l[T997l) . In general, the renormalization group flow generates spins whose average size grows without 
limit as a power law. In fact, these large spin clusters grow in a fashion which can be described as a "random walk in spin space", 
so that the total spin and the cluster size scale as 

52 - L. (75) 
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Moreover, the energy-length scale relation has a conventional power-law form n '-^ L~^. For disorder distributions which are 
not too singular, the dynamical exponent has a universal value z « 4.5 (for stronger disorder, the value of z is non-universal). 
Employing similar arguments as outlined above for the random singlet phase one finds 

X{T)-^, (76) 

Interestingly, the z exponent does not show up in the susceptibility to leading order. This is because there is a cancellation of the 
temperature dependences of the number of spin clusters and squared spin of the clusters (Eq. il5\ ) 

X (T) ^ jr^ '- ~ -. (78) 

This phase has been dubbed a "large spin phase" and has been also identified in some spin ladder systems ("Hovos and Miranda[ 
12004: Melin et al... 2002; Yusuf and Yang, 2003a b). In fact, the larger connectivity of spin ladders and higher dimensional sys- 
tems favors the formation of ferromagnetic couplings and large spins in the decimation procedure. This is easily seen if we 
consider the decimation of an antiferromagnetically coupled spin pair such that two other spins are initially antiferromagnet- 
ically coupled to the same spin. It is physically evident and easy to show that the effective coupling generated between the 
latter spins is ferromagnetic. Thus, unless the lattice geometry is fine tuned (an example of which is the two-leg spin ladder 
jHovos and Mirandal l2004 iMehn et al\ l2002t lYusuf and Vans'. '2002h\ there is a strong tendency towards a large spin phase 
formation. The singular behavior first discovered by Bhatt and Lee in higher dimensions is th us only a crossover and gives way 
to a large spin phase behavior at lower temperatures ( Lin et a/.L i2003: M otrunich et al. I l200lh . 

Generically, Griffiths phases seem to be a ubiquitous phenomenon in these random insulating magnets, usually in the vicinity 
of an infinite randomness fixed point. Systems with an Ising symmetry are especially interesting. Once again, a great deal 
of insight has been gained from an essentially exact analytical solution in the vicinity and at the quantum critical point of the 
random transverse field Ising chain due to D. S. Fisher ( Fisher, 1992. 1995.) (see also ( Shankar and M urthv. 1987)). The clean 
version of this syste m has a quantum phase transition between an Ising ferromagnet and a dis ordered p aramagnet (.Lje^^H 
ll96lVPfeutv','l97tf). The disordered version of the model is related to the McCoy-Wu model ('McCov','T969; McCov and wj, 
fl968. 1969). in which context quantum Griffiths singularities were first studied (hence the name Griffiths-McCoy singularities 
sometimes used). It can be written as 

Htfim = -J2 - E ■ (79) 

i i 

It also has a quantum phase transition tuned by the difference between the average exchange (Aj = In J) and the average 
transverse field ( A/i — In /i). Through a generalization of the Ma-Dasgupta-Hu strong disorder renormalization group procedure. 
Fisher was able to solve for the effective distributions of exchange and transverse field couplings in the quantum critical region. 
The quantum critical point itself was shown to have all the features of an IRPP, including activated dynamical scaling (with 
-0 = 1/2), diverging thermodynamic responses 

x(r)- ^ ^ ' j.'^ , (80) 



1 

Tln^Tl/T)' 



C{T)/T^ ^,_3 (81) 



and widely different average and typical correlation functions (C° {x) = {erf ■ crf^J), a — x,z) 



^"(^)-^' (82) 



lnC'^(a;) V^, (83) 

where (p = (l + v^) /2, the golden mean. Note the similarity with the result for the Heisenberg case if (p is taken to be zero. 
The singular behavior is attributed to the presence of large spin clusters at low energies which can tunnel between two different 
configurations with reversed magnetizations. The behavior off but near criticality, however, is characteristic of a Griffiths phase. 
Defining 

A. _ A , 

(84) 



var (In h) + var (In J) ' 
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where var (O) = 02 _ o is the variance of the variable O, it is found that, in the slightly disordered region (Ah ^ A j) 

xiT)^S^-^^^^, (85) 
C(T)/T^53_I_^ (86) 

whereas in the shghtly ordered region (A/j ^ Aj) 

(88) 

where z ~ 1/ (2 is a non-universal disorder-dependent expone nt, typ ical of Griffiths phases. An exact expression for z 
in terms of the bare distributions was obtained in Ref. (Igloi et all 1200 ll) . The typical and average correlations both decay 
exponentially at large distances but with different correlation lengths 

^typ ~ ^, (89) 

U-^, (90) 



such that S,av ^ S,typ and C'^ (x) 3> exp In (x). The exponent of the correlation length of the average correlation function is 
the one that satisfie s (in fact, s aturates) the bound v > 2/ d ~ 2 of Ref. dChaves et al\.\\9^^ . which is the generalization of the 
Harris criterion (.Harrislfl974l) to a disordered critical point. 

There has been independent (i.e., not based on the approximate RG scheme) numerical confirmation of Fisher's results on the 
ID random transverse field Ising model through the mapping to free fermions (Young and Rieger, 1996). Furthermore, similar 
behavior was also found in highe r dimensions by Quantu m Monte Carlo ( Pich etai, 1 998 ) and by a numerical implementation of 
the above decimation procedure jMotrunich et al. " I2OOII) . In addition, a random dilute Ising model in a tra nsverse field has been 



shown to have these same general properties close to the percolation transition in any dimension d > 1 dSenthil and Sachdevl 
[1596). 

A similar scenario h as also been discussed phenomenologically for a quantum Ising spin glass transition in higher dimensions 
(lThillandHu sei ll995h . The presence of a Griffiths phase in a model of a quantum Ising spin glass in a transverse field has been 
confirmed numerically in d = 2 and 3 (Guo et al, 1994, 1996; Rieger and Young, 1994, 1996), although the IRFP was neither 
confirmed nor ruled out at the quantum phase transition point due to the small lattice sizes used in those studies. However, 
it is believed that this quantum critical point is also governed by an IRFR The argument i s based on the IRFP nature o f the 
random transverse field Ising model quantum critical poin t in higher dimension s dMotrunic h et al. l l200lHPiche7flnil998h and 
the realization that frustration is irrelevant at the IRFP ( Motru nich et al I I200I . Indeed, at low energies the weakest link of a 
given loop of spins is almost certainly infinitely weaker than the other hnks. Therefore, one can neglect it and the loop is always 
unfrustrated. 



2. General properties of quantum Griffiths pliases 

The general nature of Griffiths phases can be elucidated in a simple and general fashion as we now show. We have seen 
in Section IT V. A . 6 1 ho w it arises in the context of the electronic Griffiths phase. In that case, the fluctuators are local moments 
coupled by a Kondo interaction Ji to the fluid of conduction electrons. The energy scale governing the local moment dynamics 
is the Kondo temperature, which depends exponentially on the combination piJi, Tki ~ D exp {—1/piJi), where pi is the local 
conduction electron density of states. We showed how the exponential tail of the distribution of the quantity l/piJi leads to a 
power law distribution of Kondo temperatures, P [Tk) ^ T^^^. This exponential tail is naturally obtained from localization 
effects. The final result is a collection of essentially independent Kondo spins whose characteristic energy scales are distributed 
in a power law fashion. 

The magnetic Griffiths phase arises in a wholly analogous fashion. In the vicinity of a quantum phase transition spatial 
fluctuations due to disorder generate droplets of the ordered phase inside a system which is on the whole in the disordered phase 
and vice-versa. Let us for simplicity focus on the disordered side of the transition. For the simplest case of uncorrected disorder, 
the statistics of these rare regions is Poissonian and the probability of a region of the ordered phase with volume L"^ is 



P {L'^) ~ exp {-cL'^) , 



(91) 
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where c is a constant tuned by the disorder strength. This droplet of ordered phase has a net magnetic moment. In the case of a 
ferromagnetic transition, the moment is proportional to the droplet size, whereas in the case of an antiferromagnet the net droplet 
moment is cx due to random incomplete spin cancellations on the surface of the droplet. These net moments are subject 

to quantum tunneling between states with reversed magnetizations. Since the driving mechanism is single-spin tunneling the 
total tunnehng rate of a droplet is exponential on the number of spins 

A-tJoexp(-6i'*) , (92) 

where b is another constant related to the microscopic tunneling mechanism and ojq is some frequency cutoff. The tunneling rate 
gives the energy splitting between the lowest and first excited states of the droplet, higher excited states lying far above in energy. 
Now, the interesting thing is that the exponentially small probability of a large cluster is "compensated" by the exponentially 
large tunneling times between states, such that the distribution of energy splittings (gaps) is given by a power law 

P (A) - [ dL'^P {L'^) 6[^-L0a exp {-bL"^)] - A"-\ (93) 



where a — c/h. We are then left with a distribution of independent droplets, whose frequencies are distributed according to 
a power law. Thus, although the microscopic origin of the magnetic and the electronic Griffiths phases are very different, the 
end result, namely, a collection of fluctuators whose energy scales are distributed in a power-law fashion, is common to both 
mechanisms. 

Given the power-law distribution of energy scales, thermodynamic and dynamical properties follow immediately. For exam- 
ple, the spin susceptibility at temperature T can be calculated by assuming that all the droplets (Kondo spins) with A > T 
{Tk > T) are frozen in the non-magnetic ground state, whereas all the others are essentially free and give a Curie response. The 
error introduced by the borderline moments is small if the distribution is broad as in Eq. ( I93> . Thus, the susceptibility is given 
by the number of free spins at temperature T, n (T) = dAP (A) times 1/T 

X(r)^^^^. (94) 

Note that, in the magnetic droplet case, we have neglected the droplet moment /i cx i*^ oc In'^^'^ A (where ip depends on the 
nature of the magnetic order), because it only gives rise to a negligible logarithmic correction to the power law (see, however, 
Section lrV.B.4> . Analogously, the total entropy S (T) cx n (T), such that the specific heat is given by 



C (T) _ dS (T) 



P(T)-^. (95) 



T dT y ' J, 

Similarly, the droplet dynamical susceptibility is 

xl<,pH~tanh(^A^5(^_A). (96) 

Averaging over the droplet distribution 

x"M^c."-itanh(^). (97) 

Interestingly, if the fluctuatorhas relaxational dynamics (as in the case of a Kondo spin) 

„ ^X(T)T{T) 

^rel (^) - ^2 + r2 (T) ' ^ ^ 

where F (T) is the relaxation rate and x [T] is the static susceptibility. Typically at high temperatures, F (T) ~ T, whereas as 
T ^ 0, F (T) - A. Since x (2^) T (T) is an almost constant smooth function of T, we can ignore it in Eq. ( I98L We thus obtain 
at low temperatures after averaging 



^'i^)-^sr^2-^"~'^ (99) 

which has the same form as Eq. (I97> . 

Similar considerations apply to the magnetization as a function of magnetic field H. The fluctuator magnetization has the 
following limiting behaviors 

H:»T,A; 

Mfi^, {H, T, A) ^ <( H/A H^T<t: A; (100) 

i7 < A < T, 
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where fi is the fluctuator momen t. This ap pHes both to a magnetic droplet in a transverse field dCastro Neto and Jonesl bOOtf) 
and to a Kondo spin (A = Tk) jHewso nl ll993..) . Averaging over A with P (A) from Eq. i93\ . it is clear that the large field 
behavior is dominated by the fluctuators with A <ti H, such that 









where Jl is an average moment (we have again neglected the possible logarithmic contribution coming from the moment size). 
At lower fields, this crosses over to a linear behavior, with a coefficient given by the (singular) susceptibility of Eq. ( I94t . 

The physical properties of Griffiths phases are therefore rather generic and are determined almost exclusively by the exponent 
of the power law distribution of energy scales, quite independent of the underlying mechanism. This makes it very difficult to dis- 
tinguish microscopic models solely on the basis of this generic power law behavior of macroscopic properties. We will, however, 
point out in Section irv.B.4l some possible measurements which can distinguish between different microscopic mechanisms. 

3. Possibility of quantum Griffitlis pliases in itinerant random magnets 

The discovery of a plethora of quantum Griffiths phases in the vicinity of phase transitions of several models of disordered 
insulating magnets immediately suggests the possibility of similar phenomena in metallic compounds. As we have seen in Sec- 
tion llll.Cl many of the systems of interest have the required ingredients, namely, the presence of extrinsic disorder and/or some 
kind of magnetic order whose critical temperature can be tuned to zero by some external parameter such as external or chemical 
pressure. Furthermore, the diverging susceptibility with tunable exponents observed close to the metal-insulator transition in 
doped semiconductors on the metallic side was found to have many simi larities to Griffiths singularities (.Bhatt and Fisher.. 19921 
iBhatt and Le^ll98llll982llDobrosavhevic et aail992HPaalanen e7aZi ri988). 

• Quantum Griffiths phase in the theory of Castro Neto, Castilla and Jones 

The possibility of a m agnetic Griffiths phase as the origin of NFL behavior in Kondo alloys was first put forward by Castro 
Neto, Castilla and Jones llCastro Neto et aZ.tll998l) . In that work, a disordered anisotropic Kondo lattice system was analyzed, 
since heavy fermion systems are known to exhibit large spin-orbit interactions which tend to strongly break Heisenberg spin 
SU(2) symmetry. Due to the competition between the RKKY interaction and the Kondo effect in a disordered environment, it was 
argued that large spatial fluctuations would give rise to regions where either of the two competing tendencies would dominate. 
The local moments that belong to a region where the local Kondo temperature is especially large are efficiently quenched at 
low temperatures and contribute to the formation of a disordered heavy Fermi liquid. On the other hand, anomalously small 
Kondo temperature fluctuations give rise to local spins which, though still subject to weak Kondo spin-flip processes, are left to 
interact with other similar partners through the RKKY interaction. Assuming strong anisotropy, the magnetic interactions lead to 
a random Ising model for the unquenched spins. Furthermore, the Kondo spin-flip terms were shown to give rise to an effective 
transverse field, i.e., if the Ising model variables are erf , anisotropic Kondo scattering generates a local field in the x-direction. 
One is then left with an effective random Ising model in a transverse field, which was discussed in Section irV B.ll As we have 
seen, there is numerical evidence that the order-disorder quantum phase transition of this model is governed by an IRFP flanked 
on both sides by Griffiths phases with non-universal tunable exponents. In the picture proposed by Castro Neto et ai, alloying 
tends to enhance quantum (Kondo) correlations, which act as a destabilizing mechanism on long-range magnetic order Thus, a 
power law distribution of energy scales is immediately obtained as is generically expected within Griffiths phases, as explained 
in Section irV.B.2l The calculation of physical quantities then follows from the arguments of that Section and they obtained, in 
terms of a non-universal tunable parameter A (which is the same as the a of Section irV.B.2> 
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(5x(T)/x(T)^r-V2^ (107) 

where x^'^-' (T) is the non-Hnear magnetic susceptibility, Xioc i'^) imaginary part of the local frequency-dependent suscep- 
tibility, Ti (T) is the NMR spin relaxation time, and Sx (T) is the root mean square susceptibility due to the disorder fluctuations. 
As we have seen, A < 1 implies NFL behavior. Some of these predictions have found experimental support in the NFL Kondo 
alloys Thi_2;U2;Pd2Al3, Yi-^jU^Pda, and VCu^-xMx (M = Pd, Pt , where the specific heat and susceptibility values could be 
well fitted to the forms in Eq. ( llQ3t with A ranging from about 0.6 to about 1 ( d e Andrade et al, 1998). Some other systems 
have been reanalyzed by Stewart in his review tiStewart„ .2001.) in light of the above results and found to be describable by such 
power laws with appropriately chosen values of A. 

Here again, the anomalous power laws are a result of a power-law distribution of energy scales for the spin fluctuators (see 
the discussion in Section ITV.B.2> . In the specific case here, in which the system teeters at the onset of magnetic order, these 
are large clusters of one phase (say, the ordered one) inside the other phase (the disordered one). Sample averaging then leads 
to the results in Eqs. (I103I107> . Although these are different from the Iow-Tr- spins of the electronic Griffiths phase, the same 
phenomenology is obtained in either case. 

• Effects of dissipation by the conduction electrons 

These early results of the magnetic Griffiths phase route towards NFL behavior did not, however, take into account the effect 
of dissipation on the tunneling of the Ising spins caused by the excitation of conduction electron particle-hole pairs of the metallic 
host. These were later incorporated in a more complete treatment by Castro Neto and Jones ( ICastro Neto and Jonesll2000() . In 
that work, the authors considered essentially the same anisotropic Kondo lattice model. The same competition between the 
RKKY interaction and the Kondo effect in a disordered environment was argued to lead, in general, to the formation of clusters 
of spins locked together by their mutual interactions and with widely distributed sizes. These clusters are able to tunnel at low 
temperatures between two different configurations with reversed signs of magnetization (in the ferromagnetic case) or staggered 
magnetization (in the antiferromagnetic case). The tunneling is caused either by residual anisotropic RKKY interactions or by 
the Kondo spin-flip scattering with the conduction electrons. The latter was called a cluster Kondo effect. The picture is thus 
similar to the previous work but now there is a distribution of cluster sizes P (iV), which is argued to be given by percolation 
theory as 

P [N) = ^ . (108) 

Here, 6* is a critical exponent and iV^ is a correlation volume (given in terms of the number of spins enclosed by it) that diverges 
as N ^ \p — Pc\ ^ close to the percolation threshold p^. The origin of percolation here is the alloying-induced growth of the 
quenched spin fluid, which acts as an inert pervasive heavy Fermi liquid. 

For a general cluster with iV spins, the hare tunneling splitting between the two configurations scales like 

Ao (TV) = c^oe-^"^, (109) 

where 7 is related to the single-spin-flip mechanism and loq is an attempt frequency. Moreover, the authors also included the 
renormalization generated by dissipation due to the conduction electrons in the cluster tunneling processes. The dissipation 
constant of a cluster was shown to scale like a (TV) = {N/Nc)'^ , where is an exponent dependent on the sp ecific ordering 
wave vecto r or tunneling mechanism. This dissipative two-level system is known to be governed by the scale jLeggett et al\. 
ll987llLesage and Saleur. 1998) 
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where A is a non-universal high-energy cutoff. The renormalization acts to suppress the tunneling rate, which should be viewed 
as an effective cluster Kondo temperature. When a (N) < 1, the cluster exhibits damped oscillations between the two con- 
figurations, whereas for a (N) > 1 dissipation completely suppresses tunneling and the cluster is frozen, behaving like a 
superparamagnetic particle. From the form of a (N) it can be seen that tunneling will cease to happen for sufficiently large 
clusters, i.e., for N > Nc, where a {Nc) — 1. The authors then p roceeded to treat this c ollection of clusters of different sizes by 
a generalization of the quantum droplet model of Thill and Huse llThill and Hus ei 1 199.5ft to the dissipative case. 

The Tk distribution for small values of Tk is a result of the exponentially rare large clusters imposed by Eq. ( I108t and the 
exponentially large tunneling times given by Eqs. ( I109> and il lOt for large values of iV < Nc- As usual (see the discussion in 
Section irv.B.2L this yields the familiar Griffiths phase power law dependence P (Tk) ^ T^^^, with a tunable A, which in turn 
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is at the origin of the singular thermodynamic properties. The new element here is the cutoff imposed by dissipation at A'^ = Nc. 
It replaces the above power law, at very low T^, by an extremely singular form 

PiTK)- ^ , ,1 . y (111) 

The power law beha vior is vahd above a crossove r scale Tk ^ T^, whereas Eq. il 1 1> is obtained for Tk < T^, where 
can be related to N^, dCastro Neto and Jonesll2000h . Thus, for scales above Griffiths effects should be observable. However, 
for temperatures below T^, more singular responses are obtained, such as x (T) ^ 1/ [Tin (A/T)] (superparamagnetism). 
The conclusion then is that the dissipation caused by a metallic environment has a dramatic effect on Griffiths phase behavior, 
effectively suppressing it at low enough temperatures and confining it to an intermediate temperature crossover regime. 

• Theory ofMillis, Morr, and Schmalian 

Another imp ortant work on the possibi lity of Griffiths singularities in metallic disordered systems is the one by Millis, Morr, 
and Schmalian llMillis ef aai200lU2002h . These authors considered the effects of disorder fluctuations on an almost critical 
system with Ising symmetry both with and without dissipation due to conduction electrons. The approach is based on a coarse- 
grained Landau action description, possibly modified by an ohmic dissipation term, in which disorder fluctuations couple to 
the quadratic term, locally chang ing the value of Tc- The authors confined themselves to situations in which the effective 
dimensionality in the Hertz sense jHertzl Il976h c?e// = d + 2 is above the upper critical dimension 4 and the analysis can be 
carried out at the mean field level (see Section lTlI.A.H . They then considered the behavior of rare disorde r fluctuation s which 
are able to nucleate a locally ordered region ("droplet") inside the paramagnetic phase. In the first work ( iMillis et al\. 1200 iT) . 
the behavior of a single droplet was analyzed. Although they considered point, line and planar defects, we will focus here on 
the case of point droplets. It was shown that the order parameter decays very slowly 1/r) in the region right outside the 
droplet core (defined approximately as the region where the local Tc (r) > T) and not further from it than the clean correlation 
length ^. This slow decay region, which can be very large close to criticality, is essential when it comes to dissipation. They 
show that these droplets will tunn el between two rever sed-magnetization configurations at a rate given by the usual expression 
for a dissipative two-level system jLeggett et adll987 t). analogous to Eq. ( II 10> . However, the dissipation constant is far greater 
than the critical value and actually diverges at criticality. Therefore, tunneling will be totally suppressed sufficiently close to 
the quantum critical point. Millis et al. clarify this result by noting that the very large effective droplet size implies that many 
conduction electron angular momentum channels will be scattered by it, as opposed to the simple s-wave case of a point scatterer. 
Similar consid erations were also made concerning the effect of ma gnetic impurities in a system close to a ferromagnetic quantum 
critical point ( Larkin and Melnikovlll97 2: Maeb ashi et al. I I200I . 

In a subsequent paper, Millis et al. extended this analysis and considered the collective effect of a distribution of droplets 
with different sizes and strengths (Millis etal, 2002). By relating the energetics and the size of a single droplet and using 
standard statistical methods, they determined the form of the droplet distribution function. Using the results for the tunneling 
rate of a single droplet the overall response of the system could be calculated. It was shown how the usual Griffiths singularities 
of an insulating (i. e., non-dissipative) nearly critical magnet (as expounded in Section IT V.B.H can be re-obtained within this 
macroscopic Landau-functional scheme. More importantly, in the metallic case, the suppression of tunneling by dissipation was 
shown to completely destroy Griffiths behavior Indeed, since dissipation effectively freezes most droplets, their contribution is 
that of essentially classical degrees of freedom, giving rise to simple superparamagnetism x {T) ^ 1/T. Close to criticality, 
the largest contribution to the thermodynamic properties comes from such frozen superparamagnetic droplets. Moreover, in the 
contribution of the unfrozen droplets (i. e., those that do tunnel) only those droplets on the verge of classicality, and consequently 
with very low tunneling rates, contribute significantly, leading again to x [T) ^ 1 /T. The results of Millis et al. thus cast some 
doubts on the relevance of a magnetic Griffiths phase picture to the physics of NFL disordered heavy fermion compounds. 

Even though differing in their detailed approach, the results obtained by Millis et al. are in general broad agreement with 
those of Castro Neto and Jones. In both cases, disorder fluctuations lead to the formation of large clusters of ordered material 
inside the disordered phase, which are able to tunnel between two magnetic configurations. Very large clusters are exponen- 
tially rare. Their tunneling rates, however, are exponentially small, leading to a power-law distribution of small energy scales 
(tunneling energy splittings) and quantum Griffiths singularities at intermediate temperatures. Finally, dissipation imposes an 
upper cutoff Nc on the fluctuating cluster sizes: above this threshold, the clusters are frozen and behave like superparamagnetic 
particles with a classical Curie response. The frozen cluster contribution eventually swamps the quantum droplets and kills the 
Griffiths singularities at the lowest temperatures. Their discrepant conclusions with regard to the applicability of their theories 
to Kondo alloys boils down to the assumptio ns made in each case about the strength of dissipation (see the recent exchange 
dCastro Neto and Jon es. 2004a b: Millis et a/.L 120 04)). Millis et al. argue that in heavy fermion systems the local moments are 
strongly coupled to the conduction electrons such that the most natural assumptions lead to a dissipation constant which is of 
the same order as the other energy scales. As a result, the largest fluctuating cluster Nc would be fairly small (of the order of a 
few sites) and the energy scale separating Griffiths behavior from superparamagnetic response would be of order of the energy 
cutoff (e.g., the Kondo temperature of the clean system). This would leave only a rather small window of temperatures in which 
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quantum Griffiths effects might be observable, effectively ruling out this mechanism as the source of NFL behavior in disordered 
heavy fermion systems. Castro Neto and Jones, on the other hand, argue that the conduction electron density of states in the 
region of the ordered cluster is not renormalized by the Kondo effect and is small. The dissipation is therefore considerably 
reduced and Nc is correspondingly enhanced, leaving a sufficiently wide range of temperatures for the Griffiths singularities 
to be detectable in heavy fermion systems. Further progress in determining which of these pictures, if any, applies to heavy 
fermion alloys will likely come from a direct determination of the strength of dissipation in real systems, either theoretically or 
from experiments. 

• Vojta theory of QCP "rounding" 

The essential gi st of the resu lts of Millis et al. can be put in a simpler, more general form and its consequences extended, as 
shown by T. Vojta ('Voita', "2003'). The argument is based on the same coarse-grained Ising-symmetry Landau approach adopted by 
Millis et al. Vojta pointed out that droplet regions at T = are bounded in the spatial directions but infinite along the imaginary 
time dimension. Since the ohmic dissipation (~ \^^n\) term of the Hertz action is equivalent to a ^ 1/ (t — r')^ kernel in the 
imaginary time direction, each droplet at T = is equivalent to one-dimensional Ising models with ^ interactions. These 
are known to undergo a phase transition at finite temperature (equivalent to a finite coup ling constant in the quantum case, which 
means a large enough droplet), i.e., they are above their lower critical dimension jThoul ess. 1969). Therefore, sufficiently large 
droplets will order (freeze) close enough to the critical point. Interestingly, this is in complete agreement with results of Millis et 
al and lends a new, simpler perspective to the rather involved calculations of that work. Considering the statistics of the droplet 
sizes for different types of disorder distributions, it was shown then that the clean quantum phase transition is "rounded" by 
disorder For a Gaussian disorder distribution for example, the magnetization in d dimensions at T = is always finite 



where g measures deviations from the clean critical point, B is a constant, and 4> is the finite-size scaUng shift exponent of 
the clean d-dimensional system. The term "rounding" here relates to the fact that the onset of order is not a collective effect 
as in conventional phase transitions but is rather the result of a sum of many finite-size frozen droplet contributions. At finite 
temperatures, droplet-droplet interactions induce a conventional thermal phase transition. If these interactions are assumed to be 
exponential in the droplet separation, the critical temperature is a double exponential of g 



where A is a constant. Finally, for non-dissipative dynamics lixin | is replaced by cj^j. The equivalent one-dimensional Ising model 
is short-ranged and thus at the lower critical dimension. This leads naturally to the Griffiths phenomena akeady identified by 
other methods. The analysis of Vojta, therefore, puts a number of different results on a common setting by means of some simple 
and transparent arguments. 

More recently, the an alysis of the Ising symmetry has been extended to the Heisenberg (in fact, continuous, O (N)) case by T. 
Vojta and J. Schmalian (Voita and SchmaliaDll200 4V Whereas the dissipative Ising droplet is above its lower critical dimension 
(when viewed as a one- dimensional system i n the time direction with 1 /r^ interactions), the Heisenberg analogue is at its 
lower critical dimension llBrund 1200 j Llovcei 1 1 9691) . Therefore, in the Vojta language, the itinerant, disordered, nearly critical 
Heisenberg magnet is analogous to the insulating, disordered, nearly critical Ising system. It is no surprise then that Griffiths 
singularities should arise in this system as well. Indeed, this is the main result of Ref. ( Voita and Schmalian, 2004), which 
shows that a power-law distribution of cluster tunneling scales is obtained in the by now familiar fashion: P (A) - A''/^-^. In 
particular, the non-universal tunable dynamical exponent z could be calculated in a controlled 1/N expansion as a function of 
disorder strength. Similar calculations in an XY magnet were also performed in Ref jLoh et aEl EoOS) with similar conclusions. 
The presence of Griffiths singularities in nearly critical metals with continuous order parameter symmetry suggests the exciting 
possibility that the associated critical point, by analogy with the insulating Ising case, is also governed by an IRFR 

• Long range RKKY interactions and the cluster glass phase 

The Griffiths phase of itinerant systems with continuous symmetry was obtained by neglecting the effe cts of drop let-droplet 
interactions. These were incorporated recently in an extended dynamical mean field theory fashion iDobrosavlievic and Mirandal 
|2005). These interactions are generically long-ranged in metallic systems due to the RKKY interactions. The presence of 
disorder is known to lead to an exponential suppression of the average RKKY interaction, because it introduces a random phase 
in the oscillating par t. However, the fluctuations around the average retain their long-ranged ~ l/r'' amplitude and are random 
in sign jjagannatha n et a/.L Il988t iNarozhnv et al.l 1200 ih . A particular droplet will therefore be coupled to a large number of 
other droplets and subject to their Weiss dynamical molecular field. A self-consistent treatment leads to a Weiss field dynamics 
set by the average local susceptibihty of all the other droplets Xioc i^n)- In the very dilute limit far from the critical point, the 
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Figur e 14 Phase diagram of an iti nerant system with disorder, close to the clean quantum critical point (QCP), following 
Ref. ( Dobrosavlievic and Mirand3 . l2005h . Shown are the clean critical line (dotted) and that in the presence of disorder if non-Ohmic dis- 
sipation is ignored (dashed). In this case, a Griffiths phase emerges close to the magnetically ordered phase. When dissipation is present, 
sufficiently large droplets "freeze", leading to the formation of the "cluster glass" phase (shaded region), preceding the uniform ordering. 

results of Vojta and Schmalian, which neglect droplet-droplet interactions, can be taken as a zeroth order approximation and an 
instability analysis can be performed. This has the usual power law form 



The Weiss dynamical field adds to the usual ohmic term \(^n\) and dominates at low frequencies if d/z < 2, which occurs 
even before the Griffiths singularity in the susceptibility (which requires d/z < 1). This sub-ohmic dissipation term places 
a single droplet above its lower critical dimension (Brutio, 2001; Joyce, 1969), leading again to freezing of sufficiently large 
droplets and eventually to ordering, in analogy with the Ising system. The conclusion is that, within this picture, long-ranged 
droplet-droplet interactions generate additional dissipation, over and above the one induced by the itinerant carriers, which 
is able to trigger magnetic order of a cluster-glass type. Interestingly, the existence of a very similar cluster glass phase was 
proposed in early work, over thirty years ago JSherrington and Mihillt 11974) . The Griffiths phase will be "hidden" inside this 
ordered phase (Fig. ll4l i. It should be stressed that this argument cannot be generalized to insulating systems, where the presence 
of short-ranged interactions casts doubts on the validity of the mean field treatment, at least in low dimensions. Likewise, for 
itinerant systems, this approach is likely to break down in low dimensions and the Griffiths phase scenario may survive. An 
interesting unsolved question is the value of the lower critical dimension for this class of disordered systems. 

• Internal quantum dynamics of droplets 

Finally, we should like to mention work which takes into account the internal quantum dynamics of droplets with a special 
geometry. Shah and Millis considered a "necklace" of spins in a ring geometry coupled as a finite XXZ chain ( Shah and Milli J 
|2003). They showed that internal quantum fluctuations of the necklace facilitate the reversal of the droplet magnetization, whose 
amplitude now decays as a power of the number of droplet spins, instead of an exponential. As a result, the nature of the droplet 
dynamics is closer to the multi-channel Kondo effect and the spin flips actually tend to proliferate at lower temperatures. The 
special geometry of a circular ring enforces a residual chira l symmetry of the necklace which protects a t \ yo-channeI fixed point, 
with its well-known anomalous properties dAffleck and Ludwig. ■1991a.b,c: .Ludwig and AfflecH Il99lt iNozieres and Blandinl 
0.980). This symmetry is not expected to be present in the case of disorder-induced droplets. Kondo-typ e dynam ics at low 
temperatures was also proposed for magnetic impurities with XY symmetry in a nearly ferromagnetic metal (iLoh et a /.L 120051) . 

4. On the applicability of the magnetic Griffiths phase phenomenology to metallic disordered systems 
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An important unsolved question that remains is whether a magnetic Griffiths phase (i.e., a Griffiths phase associated with 
a magnetic phase transition) can be realized in an observable temperature/frequency range in metallic systems. As already 
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discussed in Section IIV.B.3I a great deal of controversy surrounds the question of whether dissipation due to the conduction 
electrons is weak enough to allow for Griffiths singularities to be observable in an appreciable temperature range. However, 
a number of additional points, which are independent of the answer to this question, are also worth mentioning regarding the 
applicability of this scenario to real metallic systems. 

• Less singular thermodynamic response of heavy fermion systems: 

We have seen that power law anomalies are observed in both doped semiconductors and heavy fermion materials which 
conform, in principle, with the Griffiths phase phenomenology (Sections llll.Bl and llll.O . There are some trends, however, which 
seem to distinguish these two classes of systems. Griffiths phases are characterized by power-law divergent thermodynamic 
properties such as 



(116) 
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where the power-law exponent a is non-universal and tunable by the disorder strength. NFL behavior is signaled by a < 1. A 
significant feature of the phenomenology of the above compounds is the fact that a typically lies in the range 0.7 — 1.0 (a = 1 
impl ying a logarith mic divergence) in the case o f heavy fermion compounds (with a few cases below 0.7 but always above 0.5, 
see ( Stewart', '2001*) and ("de Andrade et flZ.','l998")), whereas doped semiconductors show a more singular response a w 0.3 — 0.4 
|Paalanen and Bhatt, 1991; Sarachik, 1995). A possible explanation to this interesting trend is the fact that localized magnetic 
moments are well formed and stable in heavy fermion systems, whereas they are induced by disorder fluctuations in doped 
semiconductors. Therefore, interactions among local moments are likely to be stronger in the former than in the latter. Such 
interactions, usually neglected in Griffiths phase theories, where the droplets are assumed dilute and independent, contribute 
to quench the local moments and may be at the origin of the less singular spin entropy available at low temperatures in heavy 
fermion materials. 

• Wilson ratio: 

As shown in Section irV.B.2l the Griffiths phase thermodynamic responses can be obtained from the scaling of a few physical 
quant ities with temperature. The results on the random quantum Ising model are very instructive in this respect llFisheilll992l 
Il995h . The spin susceptibihty, for example, is given by 

X(r) - ^ y ^ % (117) 

where /i (T) is the average value of the magnetic moment per cluster and n (T) is the number of active clusters, both at temper- 
ature T. The specific heat is obtained similarly from the entropy 

5(r) - n(T)ln2, (118) 

C{T)/T^^. (119) 

As mentioned in Section irv.B. II in the disordered Griffiths phase of the one-dimensional random Ising chain in a transverse 
field. 



It can be seen that the Wilson ratio 



^(^)^,.-..MimL, (120) 
C{T)/T^5'',^^^. (121) 

Rw^^^^-W/T)]\ (122) 

which is logarithmically divergent as T ^ 0. The origin of this divergence can be traced back to the scaling of the magnetic 
moment of the ordered (ferromagnetic) clusters with size or energy scales, which leads to /i ~ In (ilo/T). Thus, the Wilson 
ratio is a useful quantity to estimate the cluster moment at low temperatures, as the cluster number density which appears both 
in the susceptibility and the specific heat cancels out. These results are expected to be valid even in higher dimensions. Since 
we expect the cluster moments to grow with size, both in ferromagnetic and in antiferromagnetic Griffiths phases (in the latter, 
fi ^ L'^'^^^^/'^), the Wilson ratio should diverge, albeit slowly, as T ^ 0. Although perhaps hard to determine, a careful 
examination of the Wilson ratio may be a useful guide to test the applicability of the magnetic Griffiths phase scenario. 
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• Observed entropy and the size of Griffiths droplets: 

The generic quantum magnetic Griffiths phase is buih upon large clusters which tunnel between two reversed (staggered or 
uniform) magnetization states. Therefore, the available entropy per cluster is Sd — fcs In 2, the other microscopic degrees of 
freedom being effectively frozen. On the other hand, the relevant clusters are generically assumed to be large, with a typical 
number of iV^ spins per cluster If the total number of clusters is of order Nd, the total entropy can be estimated as S ^ NdSd, 
and the entropy per spin is 

S _ NdSd Sd fcBln2 



where Ngpins is the total number of spins. We conclude that the entropy per mole of spins is decreased from the typical value of 
Rln2 by a number of the order of the typical cluster size. We thus expect that conventional magnetic Griffiths behavior should 
be characterized by divergent power laws with rather small amplitudes. Conversely, if Griffiths phase behavior with a sizable 
(~ R In 2) molar spin entropy is observed, it necessarily follows that the rel evant fluctuator s invo lve a small, of order 1 , nu mber 
of spins. All the candidate heavy fermion compounds analyzed if Refs. ("Stewart', '2001) and (de Andrade et al.\ '1998*) have 
molar spin entropies of the order of R, which seems to leave little room for large clusters ( Asuiar et ai, 2003). Therefore, the 
electronic Griffiths phase (Section riV.A.6> . whose fluctuators are individual Kondo spins seems a much more viable explanation. 
Of course, if dissipation freezes clusters larger than a critical size Nc, these clusters will contribute a Curie-like term, which 
should be taken into account. However, it is clear that the above argument also limits the low temperature entropy of the frozen 
clusters, even though their contribution is more singular The conclusion that Smoiar ^ R implies small fluctuators seems 
therefore inescapable. 

• Temperature range of Griffiths phase behavior: 

Another question relates to the range of temperatures where the Griffiths singularities may be observable. We have seen that 
the strength of dissipation generated by the conduction electrons is a very important input in the determination of this range 
(see Section lTV.B.3> . Only very small dissipation rates are compatible with a wide range of Griffiths singularities. However, 
even if the dissipation is negligibly weak, the range of temperatures available for Griffiths anomalies is still fairly restricted. Let 
us consider the phenomenology of insulating disordered magnets discussed in Section lTV.B.ll As we have seen, all the known 
Griffiths phases in these systems occur in the vicinity of a phase or a point which is governed by an infinite randomness fixed 
point. Formally, this is usually seen through the divergence of the dynamical exponent z as the IRFP is approached. Now, as 
the system approaches the IRFP, there is a line of temperature crossover that approaches zero, as in any quantum critical point. 
Above the crossover line, the system is still governed by the IRFP, even off criticality. Only below the crossover line can the 
non-critical behavior, in particular the quantum Griffiths singularities, be observed ( Sachdev, 1999). The functional form of this 
crossover line can be easily obtained for systems governed by an IRFP (Fisher, 1995). The primary feature of the IRFP is the 
"activated dynamical scaling" which relates energy and length scales (see Section lTV.B.H 

17-17oexp(-L'^). (124) 

Close to criticality, there is a correlation length ^ which sets the scale beyond which Eq. ( I124> no longer holds and conventional 
dynamical scalin g is recovered. As shown by Fisher (Fisher, 1995), this is the "true" correlation length in the sense of the 
criterion of Ref. dChaves et al l ll986h : beyond this scale, the system "knows" it is non-critical and most magnetic clusters can 
be treated as effectively independent, leading to typical Griffiths singularities. This "true" correlation length diverges as a power 
law as the IRFP is approached {v 0) 

e-^, (125) 

where ly satisfies the HaiTis criterion vd > 2 ("Chaves et al. ",'198^;' Harris, '1974V This is the correlation length of the average 
correlation function. By plugging this characteristic length into Eq. ( I124> and setting — Tcross we can determine the crossover 
line 

Tcross ~ f^oexp (^-^^ • (126) 

In the particular case of the ID Ising model in a transverse field vifj = 1. We thus see that, quite generically, quantum Griffiths 
behavior is expected to occur below an energy scale that is exponentially smaller than the natural energy scales of the problem. 
Above this scale, the more singular universal behavior characteristic of the IRFP should be observed. This represents a severe 
restriction on the range of temperatures where Griffiths effects could be observed. Apparently, the diluted Ising model analyzed 
in Ref. (Senthil and Sachdev, 1996) is an exception to this kind of behavior. This is probably due to the peculiar percolating char- 
acter of the transition in that case. Since metallic systems are characterized by long-ranged RKKY interactions, it is reasonable 
to expect that percolation is less likely to play a significant role and the more generic behavior of Eq. ( I126> should apply. 
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C. Itinerant quantum glass phases and their precursors 

1 . Inherent instability of the electronic Griffiths phases to spin-glass ordering 

So far, we have discussed the electronic and the magnetic Griffiths phase scenarios for disorder-induced non-Fermi liquid 
behavior Both pictures envision the formation of rare regions with anomalously slow dynamics, which under certain conditions 
dominate the low temperature properties. Neither picture, however, seems satisfactory for the following key reason: in both 
cases the resulting NFL behavior is characterized by power law anomalies, with non-universal, rapidly varying expo nents. In 
contrast, most experimental data seem to show reasonably weak anomalies, close to marginal Fermi liquid behavior dStewartl 
I2OOD . 

Physically, it is clear what is missing from the theory. Similarly as magnetic Griffiths phases, the electronic Griffiths phase is 
characterized (Mira nda and DobrosavheviHl2001btlTanaskovic et al: . 2004) by a broad distribution P [Tk) ^ [TkT of local 
energy scales (Kondo temperatures), with the exponent a ^ rapidly decreasing with disorder W . At any given temperature, 
the local moments with Txii) < T remai n unscreened. As disorder i ncreases, the number of such unscreened spins rapidly 
proliferates. Within the existing theory ( Miranda and Dobrosavliev ic. '2001b': M iranda et a l.. 1996. U997at iTanaskovic et al\. 



12004) these unscreened spins act essentially as free local moments and provide a very large contribution to the thermodynamic 
response. In a more realistic description, however, even the Kondo-unscreened spins are not completely free, since the metallic 
host generates long-ranged Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions even between relatively distant spins. 

In a disordered metal, impurity scattering introduces random phase fluctuations in the usual periodic oscill ations of the RKKY 
intera ction, which, ho \ yeven retains its power law form (although its average value decays exponentially dJagannathan et alX 
Il98a iNarozhnv et aajgO OlM Hence, such an interaction acquires a random amplitude Jij of zero mean but finite variance 
< Jfj{R) >~ R~^'^ (|jagannathan et al, 1988; Narozhnv et al, 2001). As a result, in a disordered metallic host, a given spin 
is effectively coupled with random but long range interactions to many other spins, often leading to spin-glass freezing at the 
lowest temperatures. How this ef fect is particularly important in Griffiths phases can also be seen from the mean-field stability 
criterion (Bray and Mo or3,fl98flh for spin glass ordering, which takes the form 

Jxioc{T) = 1. (127) 

Here, J is a characteristic interaction scale for the RKKY interactions, and Xioc{T) is the disorder average of the local spin 
susceptibility. As we generally expect xioc{T) to diverge within a Griffiths phase, this arguments strongly suggests that in the 
presence of RKKY interactions such systems should have an inherent instability to finite (even if very low) temperature spin 
glass ordering. 

Similarly as other forms of magnetic order, the spin glass ordering is typically reduced by quantum fluctuations (e.g. the Kondo 
effect) which are enhanced by coupling of the local moments to itinerant electrons. Sufficiently strong quantum fluctuations can 
completely suppress spin-glass ordering even at T = 0, leading to a quantum critical point separating metallic spin glass from 
the conventional Fermi liquid ground state. As in other QCP's, one expects the precursors to magnetic ordering to emerge even 
before the transition is reached and produce non-Fermi liquid behavior within the corresponding quantum critical region. Since 
many systems where disorder-driven NFL behavior is observed are not too far from incipient spin-glass ordering, it is Ukely that 
these effects play an important role, and should be theoretically examined in detail. 

From the theoretical point of view, a number of recent works have examined the general role of quantum fluctuations in glassy 
systems, and the associated quantum critical behavior. Most of the results obtained so far have concentrated on the behavior 
within the mean-field picture (i.e., in the limit of large coordination), where a consistent description of the QCP behavior has 
been obtained for several models. In a few cases (Read et al., 1995), corrections to mean-field theory have been examined, but 
the results appear inconclusive and controversial at this time. In the following, we briefly review the most important results 
obtained within the mean-field approaches. 



2. Quantum critical behavior in insulating and metallic spin glasses 

• Ising spin glass in a transverse field 

The simplest framework to study the quantum critical behavior of spin glasses is provided by localized spin models such as 
the infinite-range Ising model in a transverse field (IMTF) with random exchange interactions Jij of zero mean and variance 
,P /N (N — > Qo is the number of lattice sites). 

HtFIM = - ^ J^J(Tl(T^ - r ^ (Tf . (128) 

ij i 

In the classical limit (F = 0), this model reduces to the well-studied Sherrington-Kirkpatrick model dMezard et a/.Ul98 ^. where 
spins freeze with random orientations below a critical temperature Tsg (r = 0) = J. Quantum fluctuations are introduced by 
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Figure 15 Generic phase diagram (following iRead ef o/lliggSh ) of the quantum critical behavior for spin glasses. The parameter r, which 
measures the quantum fluctuations, can represent the transverse field for localized spin models or the Fermi energy in metallic spin glasses. 



turning on the transverse field, which induces up-down spin flips with tunnehng rate ~ F. As F grows, the critical temperature 
7sG(r) decreases, until the quantum critical point is reached at F = Fc w 0.731J, signaling the T = transition from a 
spin-glass to a quantum-disordered paramagnetic state (Fig.ll5>. 

Similarly as in DMFT theories for electronic systems, such infinite range models can be formally reduced to a self-consistent 
solution of an app ropriate quantum impurity problem, as first discussed in the context of quantum spin glasses b y Bray and Moore 
( iBrav and Moorel [19801) . Early work quickly established the phase diagram jPobrosavlievic and Strattlll987h of this model, but 
the dynamics near the quantum critical point proved more difficult to unravel, even when the critical point is approached form 
the quantum-disordered side. Here, the problem reduces to solving for the dynamics of a single Ising spin in a transverse field, 
described by an effective Hamiltonian jMiller and Huselll993h of the form 

H^-^J'^ J J drdTV^(T)x(r-r')CT'(r')+F J dra^r). 

Physically, the interaction of the considered spin with the spin fluctuations of its environment generates the retarded interaction 
described by the "memory kernel" xi''' ~ ^')- appropriate self-consistency condition relates the memory kernel to the 
disorder-averaged local dynamical susceptibility of the quantum spin 



x(t-t') = (ra-(rV(T')). 

A complet e solution of t he quantum critical behavior can be obtained, as first established in a pioneering work by Miller and Huse 
(tviifler an dHuselfT993l) . These authors have set up a diagrammatic perturbation theory for the dynamic susceptibility, showing 
that the leading loop approximation already captures the exact quantum critical behavior, as the higher order corrections provide 
only quantitative renormalizations. The dynamical susceptibility takes the general form 

xK) = Xo + K + a2)'/', 

where the local static susceptibility Xo remains finite throughout the critical regime, and the spin excitations exist above a gap 

A - (r/|lnr|)^/^ 

which vanishes as the transition is approached from the paramagnetic side (here r = (F — Fc) /Fc measures the distance from 
the critical point). 

• The quantum spin-glass phase and the replicon mode 



The validity of this solution was confirmed by a generalization ( lYe et aZ.|,ll993l) to the AZ-component rotor model (the Ising 
model belongs to the same universality class as the M ~ 2 rotor model), which can be solved in closed form in the large M 
limit. This result, which proves to be exact to all orders in the 1 /M expansion, could be extended even to the spin glass phase, 
where a full replica symmetric solution was obtained. Most remarkably, the spin excitation spectrum remains gapless (A — 0) 
throughout the ordered phase. Such gapless excitations commonly occur in ordered states with broken continuous symmetry, but 
are generally not expected in classical or quantum models with a discrete symmetry of the order parameter In glassy phases (at 
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least within mean-field solutions), ho wever, gapless e xcitations generically arise for both classical and quantum models. Here, 
they reflect the marginal stability tMezard et al\. 11986.) found in the presence of replica symmetry breaking, a phenomenon 
which reflects the high degree of frustration in these systems. The role of the Goldstone mode in this case is played by the 
so-called "replicon" mode, which describes the collective low energy excitations characterizing the glassy state. 

A proper treatment of the low energy excitations in t his regime requ ires special attention to the role of replica symmetry 
breaking (RSB) in the T — )■ limit. The original work suggeste d that RSB is suppres sed at T = 0, so that 

the simpler replica symmetric solution can be used at low temperatures. Later work ("George s et a/.ll200ll) . however, established 
that the full RSB solution must be considered before taking the T — > limit, and only then can the correct form of the leading 
low temperature corrections (e.g., the linear T-dependence of the specific heat) be obtained. 

• Physical content of the mean-field solution 

In appropriate path-integral language llDobrosavlievic and Strattlll987tf^ ef aZ.L fl993l) . the problem can be shown to reduce 
to solving a one dimensional classical Ising model with long-range interactions, the form of which must be self-consistently 
determined. Such classical spin chains with long range interactions in general can be highly nontrivial. Some important examples 
are the Kondo problem (Anderson and Yuval, 1969; Anderson et ai, 1970; Yuval and Anderson, 1970), and the dissipative two- 
level system (Leggettef a/. , 1987), both of which map to an Ising chain with l/r^ interac tions. Q uantum phase transitions 
in these problems correspond to the Kosterlitz-Thouless transition found in the Ising chain dKosterlitz. .1976) . the description 
of which required a sophisticated renormalization-group analysis. Why then is t he solution of th e quantum Ising spin glass 
model so simple? The answer was provided in the paper by Ye, Sachdev, and Read dYe et a/. ('l993). which emphasized that the 
critical state (and the RSB spin-glass state) does not correspond to the critical point, but rather to the high-temperature phase of 
the equivalent Ising chain, where a perturbative solution is sufficient. In Kondo language, this state corresponds to the Fermi- 
liquid solution characterized by a finite Kondo scale, as demonstrated by a quantum Monte-Carlo calculation of Rozenberg and 
Grempel (Rozenberg and Grempel, 1998), which also confirmed other predictions of the analytical theory. 

From a more general perspective, the possibility of obtaining a simple analytical solution for quantum critical dynamics has a 
simple origin. It follows from the fact that all corrections to Gaussian (i.e. Landau) theory are irrelevant above the upper critical 
dimension, as first established by the Hertz-Millis theory (Hertz, 1976; Millis, 1993) for conventional quantum criticalit y. The 
mean- field models become exact in the limit of infinite dimensions, hence the Gaussian solution of Refs. dMiller and HusJfl993L 
I Ye et ai. .1993) becomes exact. Leading corrections to mean-field theory for rotor models were examined by an e-expansion 
below the upper critical dimension dc = 8 for the rotor models by Read, Sachdev, and Ye, but these studies found run-away 
flows, presumably indicating non-perturbative effects that require more sophisticated theoretical tools. Most likely these include 
Griffiths phase phenomena controlled by the infinite randomness fixed point, as already discussed in Section lTV.B.ll 

• Metallic spin glasses 

A particularly interesting role of the low-lying excitations associated with the spin-glass phase is found in metallic spin glasses. 
Here the quantum fluctuations are provided by the Kondo coupling between the conduction electrons and local moments, and 
therefore can be tuned by controlling the Fermi energy in the syste m. The situation is again th e simplest for Ising spins where 
an itinerant version of the rotor model of Sengupta and Georges (Sensupta and Georges, '1995^ can be considered, and similar 
results have obtained for the "spin-density glass" model of Sachdev, Read, and Oppermann ( Sachdev et ai, 1995). The essential 
new feature in these models is the presence of itinerant electrons which, as in the Hertz-Millis approach ( HertJ 1 19761 iMillisl 
Il993h . have to be formally integrated out before an effective order-parameter theory can be obtained. This is justified provided 
that the quasiparticles remain well defined at the quantum critical point, i.e. that the quasiparticle weight Z ~ Tr- remains finite 
and the Kondo effect remains operative. The validity of these assumptions is by no means obvious, and led to co nsiderable 
controversy before a detailed quantum Monte Carlo solution of the model became available dRozenberg and Grempei ,199% . 
confirming the proposed scenario. 

Under these assumptions, the theory can again be solved in closed form, and we only quote the principal results. Physically, 
the essential modification is that the presence of itinerant electron now induces Landau damping, which creates dissipation for 
the collective mode. As a result, the dynamics is modified, and the local dynamic susceptibility now takes the following form 

XK) =Xo + (kn|+^*)'^'. (129) 

The dynamics is characterized by the crossover scale ~ r (r measures the distance to the transition) which defines a crossover 
temperature T* ^ uj* separating the Fermi liquid regime (at T ^ T*) from the quantum critical regime (at T ^ T*). At the 

1/2 

critical point xi^n) — Xo + \<^n\ , leading to non-Fermi liquid behavior of all physical quantities, which acquire a leading 
low-temperature correction of the T^f^ form. This is a rather mild violation of Fermi liquid theory, since both the static spin 
susceptibility and the specific heat coefficient remain finite at the QCR A more interesting feature, which is specific to glassy 
systems, is the persistence of such quantum critical NFL behavior throughout the metalhc glass phase, reflecting the role of the 
replicon mode. 
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3. Spin-liquid beliavior, destruction of the Kendo effect by bosonic dissipation, and fractionalization 

• Quantum Heisenberg spin glass and the spin-liquid solution 

Quantum spi n glass behavior proves to be much more interesting in the case of Heisenberg spins, where the Berry phase term 
plays a highly nontrivial role, completely changing the dynamics even withi n the paramagnetic ph ase. While the 
existence of a finite temperature spin-glass transition was established even in early work ( Brav and Moo rel ll98CII) . solving for 
the details of the dynamics proved difficult until the remarkable work of Sachdev and Ye (Sachdev and Ya. ll993l) . By a clever 
use of large- methods, these authors identified a striking spin-liquid solution within the paramagnetic phase. In contrast to the 
nonsingular behavior of Ising or rotor quantum spin glasses, the dynamical susceptibility now displays a logarithmic singularity 
at low frequency. On the real axis it takes the form 

7r 

X{uj) - ln(l/ \uj\) + i-^sign{uj). 

A notable feature of this solution is that it is precisely of the form postulated for "marginal Fermi liquid" phenomenology 
( Varma et al., 1989) of doped cuprates. The specific heat is also found to assume a singular form C ~ VT, which was shown 
I Georges et al. , 2001 ) to reflect a nonzero extensive entropy if the spin liquid solution is extrapolated to T = 0. Of course, the 
spin-liquid solution becomes unstable at a finite ordering temperature, and the broken symmetry state has to be examined to 
discuss the low temp erature propert ies of the model. 

Subsequent work ( iGeorges et ali 1200 ll) demonstrated that this mean-field solution remains valid for all finite N and gener- 
alized the solution to the spin-glass (ordered) phase. A closed set of equations describing the low temperature thermodynamics 
in the spin glass phase was obtained, which was very recently re-examined in detail (Camiavi and Rozenberg, 2003) to reveal 
fairly complicated behavior 

• Metallic Heisenberg spin glasses and fractionalization 

Even more interesting is the fate of this spin Uquid solution in itinerant systems, where an additional Kondo coupling is added 
between the local moments and the conduction electrons. The mean-field approach can be extended to this interesting situation 
by examining a Kondo-Heisenberg spin glass model (Burdin et al., 2002) with the Hamiltonian 

Hkh = -tYl + '^■) + ^if ^ S, • s, + ^ S, • S, . (130) 

(ij)CT i (ij) 

1/2 

In the regime where the scale of the RKKY interaction J = { Jfj ) is small compared to the Kondo coupling Jk, one expects 
Kondo screening to result in standard Fermi liquid behavior In the opposite limit, however, the spin fluctuations associated with 
the retarded RKKY interactions may be able to adversely affect the Kondo screening, and novel metallic behavior could emerge. 
This intriguing possibility can be precisely investigated in the mean-field (infinite range) limit, where the problem reduces to a 
single-impurity action of the form (Burdin et al., 2002) 

^i// -^T. t (r) {dr~P + V,) (r) 

dT 
Jo 



t'J2 I dT I dT'cl{T)G,{T-T')c^{T') 
a 

r0 



+ Jk I dTS{T) ■ s(r) 

"'0 



— / dr / dT'x(T-r')S(r).S(T'). (131) 
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Such a sin gle-impurity action (1131^ describes the so-called Bose-Fermi Kondo (BFK) impurity model dSenguptal l2000t 
ISi and Smifli Il996t IZarand and Demled l2002t IZhu and Sii Eool where, in addition to the coupling to the fermionic bath of 
conduction electrons, the Kondo spin also interacts with a bosonic bath of spin fluctuations, with local spectral density x (^n)- 
Because the sa me BFK m odel also appears in "extended" DMFT theories (Smith and Si, 2000) of q uantum criticality in clean 
systems dSengu pta. 2000t lSi et g/.tboOlHsTand Smiflj. [1993 IZarand and Demleil l2002HZhu and^it l2002i) . its properties have 
been studied in detail and are by now well understood. 

In the absence of the RKKY coupling (J = 0), the ground state of the impurity is a Kondo singlet for any value of Jk 7^ 0. 
By contrast, when J > 0, the dissipation induced by the bosonic bath tends to destabilize the Kondo effect. For a bosonic bath 
of "Ohmic" form (x (i^n) = Xo — C |ijJ„|), this effect only leads to a finite decrease of the Kondo temperature, but the Fermi 
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T 




RKKY coupling J 

Figure 16 Phase diagram of the Bose-Fermi Kondo model in the presence of a sub-Ohmic bosonic bath ^SenguDtaL l200(l ISi and Smithl 
[l996 ; Zarand and Demler, 2002 ; Zhu and Si. 2002). Kondo screening is destroyed for sufficiently large dissipation (RKKY coupling to spin 
fluctuations). 

liquid behavior persists. In contrast, for "sub-Ohmic" dissipation (x (wn) = Xo — C |w„|^~'^ with e > 0) two different phases 
exist, and for sufficiently large RKKY coupling the Kondo effect is destroyed. The two regimes are separated by a quantum 
phase transition (see Fig.ll6>. 

Of course, in the considered Kondo lattice model with additional RKKY interactions, the form of the bosonic bath x {^n) is 
self-consistently determined, and can take differe nt forms as the RKK Y coupling J is increased. The model was analytically 
solved within a large N approach by Burdin et al. teurdi n et a/.Ll2003) . who calculated the evolution of the Fermi liquid coher- 
ence scale T*, and the corresponding quasiparticle weight Z in the presence of RKKY interactions. Within the paramagnetic 
phase both T*{J) and Z{J) are found to decrease with J until the Kondo effect (and thus the Fermi hquid) is destroyed at 
J — fa lOT^ (here is the J = Kondo temperature), where both scales vanish ( iBurdin et all 12002: iTanaskovic et all 
|2005). At T > T*{J) (and of course at any temperature for J > Jc) the spins effectively decouple from conduction electrons, 
and spin liquid behavior, essentially identical to that of the insulating model, is established. Thus, sufficiently strong and frus- 
trating RKKY interactions are able to suppress Fermi liquid behavior, and marginal Fermi liquid behavior emerges in a metallic 
system. This phenomenon, corre sponding to spin-charge separation resulting from the destruction of the Kondo effect, is some - 
times called "fractionalization" dColeman and An drei. 1987: D emler et al !. 2002; Kagan et al, 1992; Senthil et al., 2003, 20041). 
Such behavior has often been advocated as an appealing scenario for exotic phases of strongly correlated electrons, but with 
the exception of the described model, the re are very few well establish ed results and model calculations to support its validity. 
Finally, we should mention related work dParcoUet and Georgeslll99'^ on doped Mott insulators with random exchanges, with 
many similarities with the above picture. 

We should note, however, that this exotic solution is valid only within the paramagnetic phase, which is generally expected 
to become unstable to magnetic (spin glass) ordering at sufficiently low temperatures. Since fractionalization emerges only for 
sufficiently large RKKY coupling (in the large N model Jc ~ IQT^), while in general one expects magnet;ic ordering to take 
place already at J ~ (according to the famous Doniach criterion jPoniachl Il977t) ). one expects ( iBurdin et all 12002 ) the 
system to magnetically order much before the Kondo temperature vanishes. If this is true, then one expects the quantum critical 
behavior to be very similar to metallic Ising spin glasses, i.e., to assume the conventional Hertz-Millis form, at least the for 
mean-field spin glass models we discussed here. The precise relevance of this paramagnetic spin liquid solution thus remains 
unclear, at least for systems with weak or no disorder in the conduction band. 

• Fractionalized two-fluid behavior of electronic Griffiths phases 

The situation seems more promising in the presence of sufficient amounts of disorder, where the electronic Griffiths phase 
forms. As we have seen in Sec. IIV.A.6I here the disorder generates a very broad distribution of local Kondo temperatures, 
making the system much more sensitive t o RKKY interactions. This mechanism has recently been studied within an extended 
DMFT approach ( T anaskovic et a4l200 5*). which is able to incorporate both the formation of the Griffiths phase, and the effects 
of frustrating RKKY interactions leading to spin-glass dynamics. At the local impurity level, the problem still reduced to the 
Bose-Fermi Kondo model, but the presence of conduction electron disorder qualitatively modifies the self-consistency conditions 
determining the form of x {'-^n)- 
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To obtain a sufficient condition for decoupling, we examine the stability of the Fermi liquid solution, by consider- 
ing the limit of infinitesimal RKKY interactions. To leadi ng order we repl ace xi^n) — > Xo i^n) = xi'^n', J = 0), 
and the calculation reduces to the "bare model" of Ref. (Tanaskovic et aZ.L l2Q05h . In this case, from Eq. ( 16 1> . P {Tk) ~ 7/f^> 
where a ^ and 

XoM- j dTKP{TK)x{^n,TK)^Xo{Q)-CM^~\ (132) 

where e = 2 — a. Thus, for sufficiently strong disorder (i.e. within the electronic Griffiths phase), even the "bare" bosonic bath 
is sufficiently singular to generate decoupling. The critical value of W will be modified by self-consistency, but it is clear that 
decoupling will occur for sufficiently large disorder. 

Once decoupling is present, the system is best viewed as composed of two fluids, one made up of a fraction n of decoupled 
spins, and the other of a fraction (1 — n) of Kondo screened spins. The self-consistent x (^n) acquires contributions from both 
fluids 

X = nxdc (w„) + (1 - n) Xs (^«) • (133) 
A careful analysis (iTanaskovic et a/.Ll2005l) shows that, for a bath characterized by an exponent e 

Xdc K) ~ Xdc (0) - C |w„| ; (134) 
Xs K) ^ Xs (0) - C |a;„|^-(2-^-^/'^) , (135) 

where v = v {e)\s a. critical exponent governing how the Kondo scale vanishes at the quantum critical point of the Bose-Fermi 
model. Since v > {), the contribution of the decoupled fluid is more singular and dominates aMower frequencies. Self- 
consistency then yields e = 1, as in the familiar spin liquid state of Sachdev and Ye dSachdev and Yd. Il993h . For e = 1, the 
local susceptibility is logarithmically divergent (both in ujn and T). This does not necessarily mean that the bulk susceptibility, 
which is the experimentally relevant quantity, behaves in the same manner (Parcollet and Georaes, 1999). More work remains 
to be done to determine the precise low temperature form of this and other physical quantities, and to assert the relevance of this 
mechanism for specific materials. 

As in the case where conduction electron disorder is absent, the spin liquid state is unstab le towards spin-^ l ass or dering at 
sufficiently low temperatures. However, numerical estimates for the Griffiths phase model (Tanaskovic etal., 2005) suggest 
a surprisingly wide temperature window where the marginal behavior should persist above the ordering temperature. Fig. 1171 
represents the predicted phase diagram of this model. For weak disorder the system is in the Fermi liquid phase, while for 
W > Wc the marginal Fermi liquid phase emerges. The crossover temperature (dashed line) delimiting this regime can be 
estimated from the frequency up to which the logarithmic behavior of the local dynamical susceptibility x(*^) is observed. 
The spin glass phase, obtained from Eq. ( I127> . appe ars only at t he low est temperatures, well below the marginal Fermi liquid 
boundary. Interestingly, recent experiments (MacLaughlin et a/.Ll200ll) have indeed found evidence of dynamical spin freezing 
in the milliKelvin temperature range for the same Kondo alloys that display normal phase NFL behavior in a much broader 
temperature window. 

The two-fluid phenomenology of the disordered Kondo lattice we have described above is very reminiscent of earlier work 
on the clean Kon do lattice, where the conduction electrons effectively decouple from the local moments, the latter forming a 
spin liquid state ('Coleman and Andrei', '1987"; 'Demler et al.', '2002'; 'Kagan et al.', '1992'; 'Sen thil et ali 12003', "2004). The major 
difference between the results presented in this Section and these other cases is that here local spatial disorder fluctuations lead 
to an inhomogeneous coexistence of the two fluids, as each site decouples or not from the conduction electrons depending on its 
local properties. The discussed mean-field models should be considered as merely the first examples of this fascinating physics. 
The specific features of the spin liquid behavior that was obtained from these models may very well prove to be too restrictive 
and perhaps even inaccurate. For example, the specific heat enhancements may well be overestimated, reflecting the residual 
T = entropy of the mean-field models. Nevertheless, the physics of Kondo screening being destroyed by the interplay of 
disorder and RKKY interactions will almost certainly play a central role in determining the properties of many NFL systems, 
and clearly needs to be investigated in more detail in the future. 

4. Electron glasses, freezing in the charge sector, and the quantum AT line 

Glassy behavior in disordered electronic systems is not limited to phenomena associated with the freezing of spins. In fact, 
glassy physics in the charge sector was already envisioned a long time ago in the pioneering works of Efros and Shklovskii 
fEfros and Shklovskii, 1975; Shklovskii and Efros, 1984) and Pollak (Pollak, 1984) on disordered insulators. It is expected to 
result from the competition of the long-ranged Coulomb interaction and disorder The physical picture that has emerged from 
these works is easy to understand. In absence of disorder, the Coulomb repulsion tends to keep the electrons as far from each 



47 



0.02 r 



0.015- 



H 0.01- 




0.005 ■ 



n 

05 BARE GRIFFITHS 
PHASE 



^ 0.2 



FERMI LIQUID ' 



MARGINAL FERMI 
LIQUID 



0.2 



w 



"SPIN GLASS 

0.4 



Figure 17 Phase diagram of the electronic Griffiths phase model with RKKY interactions iTanaskovic et a.l\ l200.'j) . The inset shows the 
fraction of decoupled spins as a function of the disorder strength W . 
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Figure 18 Phase diagram of the three-dimensional classical Coulomb glass, as obtained from the EDMFT approach, following Ref. 
iPankov md Dobrosavlievic. 2005). The full horizontal line indicates the glass transition temperature, and the dotted line shows where the 
entropy of the fluid solution turns negative. The inset shows the temperature dependence of the screening length as a function of temperature. 



other as possible. If quantum or thermal fluctuations are sufficiently small, the electrons tend to assume a periodic pattern, 
leading to the formation of a charge density wave (e.g. a Wigner crystal). In contrast, disorder tends to arrange the electrons 
in a random fashion, opposing the periodic pattern favored by the Coulomb interaction. When both effects are comparable the 
system is frustrated: there now exist many different low energy configurations separated by potential barriers, leading to glassy 
dynamics. 

Of course, these glassy effects are most pronounced in disordered insulators, where quantum fluctuations are minimized by 
electron localization. Over the last thirty years a large number of theoretical studies have concentrated on the physics of the 
Coulomb glass. In this review we will not attempt to describe in detail this large body of work, since many of these results 
do not directly touch upon the non-Fermi liquid physics observed in metals. We will limit our attention to those works that 
concentrated on the manifestati ons of el ectronic glass behavior on the metallic side of the metal-insul ator transition . The original 
works of Efros and Shklovskii llEfros an d Shklo vskiiIfl97llShklovskii and EfrosHl984 and PoUak llPollaklll984i) concentrated 
on classical models for the electron glass, as did the work of many followers. Because most workers used numerical approaches 
to investigate the problem, it proved difficult to incorporate the effects of quantum fluctuations into this picture, which cannot be 
avoided close to the metal-insulator transition. 

More recently, extended DMFT approaches (EDMFT) llChitra and Kothaill200(llSmith an d Sil '20001) were shown to capture 
many relevant aspect of the classical Coulomb glasses, such as the formation of the Coulomb gap ( Pankov and Dobrosavlievid 
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Figure 19 Global phase diagram for the disordered Hubbard model iPobrosavlievic ef fl/.Ll2003h . as a function of the hopping element t and 
the disordered strength W , both expressed in units of the on-site interaction U . The size of the metallic glass phase is determined by the 
strength of the inter-site Coulomb interaction V . 



l20Q5t IPas|:9r ^n(^ D9lpr9^^Yliwi d|ip^^ (Fig.ll8>. In addition, th ese theories were able to discuss the effects of quantum i 
tuations (iDobrosavhevic et al !. 2003; Pastor and Dobrosavlievi c, 1999) in the vicinity to the metal-insulator transition. Ir 



fluc- 
In the 

following we briefly describe the physical picture of the Coulomb glass emerging from these theories. 



• Coulomb gap, the replicon mode, and self-organized criticality 

The classical Coulomb glass model of Efros and Shklovskii (ES) llEfros and ShklovskiiIll975llShklovskii and Efrosill984l) is 
given by the Hamiltonian 

1 

if = ^«,n, + -^— (n, -X)(nj-if). (136) 



Here = 0, 1 is the electron occupation number, and Vi is a Gaussian distributed random potential of variance W'^ . In their 
classic work on this model ES presented convincing evidence that at T = a soft "Coulomb gap" emerges in the single-particle 
density of states (DOS) which, in arbitrary dimension d, takes the form 



g{e) ^ e""-'. (137) 

>From a general point of view this result is quite surprising. It indicates a power-law distribution of excitation energies, i.e. 
the absence of a characteristic energy scale for excitations above the ground state. Such behavior is common in models with 
broken continuous symmetry, where it reflects the corresponding Goldstone modes, but is generally not expected in discrete 
symmetry models, such as the one used by ES . Within the EDM FT theory of the Coulomb glass llMuller and loffel l2004t 
iPankov and Dobrosavlievic, 2005; Pastor and D obrosavlievid fl999h this puzzling feature finds a natural explanation: it reflects 
the emergence of the soft "replicon" mode, similarly to the classical and quantum spin-glass models we have already discussed. 
As a result, again in dir ect analogy with spin glass models, the order ed state of the Cou lomb glass is expected to display 
self-organized criticality jPazmandi et a/.LIl999h and hysteresis behavior jPastor et a/.Ll2002l) characterized by avalanches on all 
scales. 

Another interesting feature of the Coulomb glass is w orth mentioning. A s explicitly demonstrated by examining the broken 
replica symmetric solution of the mean-field equations llMuUer and loff e. 2004), the screening length is found to diverge as 
Iscr ^ 1/T as a result of the vanishing zero-field cooled compressibility ( Pa stor and Dobrosavlievid 1 199 9) within the glassy 
phase. This result is significant because it explains the absence of screening in the ground state of the Coulomb glass, in 
agreement with the assumptions of the ES theory. It also indicates the divergence of the effective coordination number in the 
T — > limit, giving further credence to the predictions of the mean-field approach. 



• Quantum Coulomb glass and the metallic glass phase 
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When quantum fluctuations are introduced in the model by including hopping amplitudes between lattice sites , one 
must examine the T = tran sition where g lassy freezing of electrons is suppressed. Such quantum melting of the electron 
glass was examined in Refs . llDobrosavliev ic ef a/. . 2003: Pastor and Dobrosavlievic. 1999). Interestingly, the theory shows 
(iDobrosa vlievic et a l.'. 2003) that the approach to an Anderson-like insulator is a singular perturbation to the stability of the 
glassy phase. Physically, this reflects the fact that the quantum fluctuations in this model arise from the mobility of itinerant 
electrons - which is expected to vanish precisely at the metal-insulator transition. As a result, the glass transition is predicted to 
generically precede the metal-insulator transition (Fig.ll9>, giving rise to an intermediate metallic glass phase, consistent with 
some recent experiments (Bogdanovich and Popo vie, 2002,). 

• Critical behavior along the quantum AT line 

As we can see from Eq. jl36> . the Coulomb glass Hamiltonian is essentially identical to that of an antiferromagnetic random 
field Ising model (RFIM) with long range interactions. The random potential Vi plays the role of a random symmetry-breaking 
field, and as a result, the finite temperature glass transition in our model assumes the character of a de-Almeida-Thouless (AT) 
line dde Almeida and Thoulessl fl975 iMezard eTall Il986h . The T = transition separating the Fermi liquid from a NFL 
metallic glass phase thus assumes the character of a "quantum AT line". A complete descripti on of this quantum critical point 
can b e obtained within the mean-field formulation, which has been studied in detail in Ref. jPahdovich and Dobrosavlievid 
l2002h . Similarly as in metallic spin glasses (see Section irV.C.2t . within the quantum critical regime and the entire metallic glass 
phase, the resistivity acquires a non-Fermi liquid T^/"^ form, as seen in some experiments (Bogdanovich and PoDovic, 2002). In 
contrast to the metallic spin glass scenario, the approach to the quantum AT line is characterized ( IDahdovich and Dobrosavlievl3 
12002) bv a crossover scale 

T* - (138) 

where St = {t — tc) / tc measures the distance to the critical point. This indicates a much broader quantum critical region in this 
case, which may explain the apparent sup pression of weak localiz ation (Fermi liquid) corrections in two dimensional electron 
systems near the metal-insulator transition ([Abrahams et al.i\2Q()ll . 



V. CONCLUSIONS AND OPEN PROBLEMS 

We have seen throughout this review that many routes can lead to disorder-induced NFL behavior. In this concluding section, 
we would like to pause to analyze general trends of what has been done so far, what is still missing, and what are the most 
promising directions for future work. 

We have seen how the phenomenology of quantum Griffiths singularities seems to afford a natural description for many of 
the anomalies seen in several systems, from doped semiconductors to disordered heavy fermion systems. On the other hand, we 
have listed several outstanding points which still need to (and should) be clarified. More generally, however, we would like to 
stress a different question. Griffiths phases rely on the existence of rare disorder configurations which generate essentially non- 
interacting fluctuators. It is the slow quantum dynamics of these fluctuators which eventually leads to the known singularities. 
This situation is often described as being one in which quantum mechanics (tunneling) is a "dangerously irrelevant operator" 
responsible for changing the scaling dimension of otherwise classical (and non-interacting) objects. On the other hand, we 
have described how droplet-droplet interactions can dramatically change this picture and be one possible destabilizing factor of 
Griffiths physics (Section [lV.B.3> . This effect represents essentially the delocalization of the local (droplet) spin modes which 
are at the heart of the Griffiths phases. An important outstanding issue in this context is the determination of the lower critical 
dimension for spin mode localization. This is especially important in the case of metallic systems where the spin interactions 
are long-ranged. 

At or below this lower critical dimension, there appears the question of the existence of an infinite randomness fixed point in 
metallic systems with continuous symmetry. These represent the best candidates so far for a quantum Griffiths phase which can 
survive down to zero temperature, as discussed in Section irV.B.3l Below the (still unknown) lower critical dimension for spin 
mode localization, and based on the well understood case of insulating systems, we would expect that this Griffiths phase would 
be a precursor of a phase transition governed by an infinite randomness fixed point. In fact, this is probably the only known route 
for such a point that does not require fine tuning in systems with continuous symmetry. Indeed, as explained in Section lTV.B.ll 
except for a few special model systems (the nearest-neighbor random antiferromagnetic Heisenberg and XX chains being the 
most prominent examples), almost all other perturbations, such as increased connectivity and/or ferromagnetic interactions, seem 
to lead to a large spin phase. These include the important case studied by Bhatt and Lee in the context of doped semiconductors. 
Therefore, the rigorous establishment of the stability of a Griffiths phase and an associated infinite randomness fixed point in 
systems with continuous symmetry and dissipation would be an important achievement. 

Above the lower critical dimension, we expect glassy dynamics to take center stage. Here, very interesting possibilities exist 
for non-trivial physics even in the paramagnetic phase, as hinted by the results and references quoted in Section irV.C.3l Other 
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questions which arise are: (i) Is there a Griffiths phase close to the spin glass quantum critical point? (ii) Is it governed by an 
infinite randomness fixed point or is it characterized by finite disorder? (iii) How is the charge transport affected by the spin 
dynamics? (iv) What is the role of dissipation? 

The exploration of glassiness in the charge sector also deserves further attention. Here again long range Coulomb interactions 
play a crucial role. The interplay with the spin sector is still largely unexplored. This seems especially important in the case of 
the disorder-induced metal-insulator transition both in 3 and in 2 dimensions. 

We hope that addressing and answering some of these questions will provide important clues as to the origin of non-Fermi 
liquid behavior in disordered systems. 
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